PM/00-03 
February 2000 



THE HIGGS WORKING GROUP: 
Summary Report 

Conveners: 

A. Djouadi\ R. Kinnunen^, E. Richter-Was^''' and H.U. Martyn^ 

Working Group: 

K.A. Assamagan*^, C. Balazs^, G. Belanger^, E. Boos^, F. Boudjema^, 
M. Drees^°, N. Ghodbane^\ M. Guchait^, S. Heinemeyer^, V. Ilyin^, 
J. Kalinowski^^, J.L. Kneur^ R. Lafaye^, D.J. Miller^, S. Moretti^^, 

M. MUHLLEITNER^ A. NlKITENKO^'^, K. OdAGIRI^^, D.P. ROY^^, M. SpIRA^^ 

K. Sridhar^^ and D. Zeppenfeld^'''^'^. 



1 LPMT, Universite Montpellier II, F-34095 Montpellier Cedex 5, France. 
^ Helsinki Institute of Physics, Helsinki, Finland. 
3 CERN, IT Division, 1211 Geneva 23, Switzerland. 
^ Institute of Computer Science, Jagellonian University, and Institute of Nuclear Physics, 

30-059 Krakow, ul. Nawojki 26a, Poland. 
5 DESY, Notkestrasse 85, D-22603 Hamburg, Germany 
^ Hampton University, Hampton, VA 23668, USA. 
^ Department of Physics and Astronomy, University of Hawaii, Honolulu, HI 96822. 
s LAPP, BP 110, F-74941 Annecy le Vieux Cedex, France. 
^ Institute of Nuclear Physics, MSU, 11 9899 Moscow, Russia. 
Physik Department, TU Miinchen, James Franck Str., D-85748 Garching, Germany. 

IPNL, Univ. Claude Bernard, F-69622 Villeurbanne Cedex, France. 
^■^ Institute of Theoretical Physics, Warsaw University, PL-00681 Warsaw, Poland. 
13 Rutherford Appleton Laboratory, Chilton, Didcot, Oxon 0X11 OQX, U.K. 
1^ Theoretical Physics Department, TIFR, Homi Bhabha Road, Bombay 400 005, India. 
1^ II. Inst. Theor. Physik, Universitat Hamburg, D-22761 Hamburg, Germany. 

1^ CERN, Theory Division, CH-1211, Geneva, Switzerland. 
1'' Department of Physics, University of Wisconsin, Madison, WI 53706, USA. 



Report of the HIGGS working group for the Workshop 
'Physics at TeV Colliders" , Les Houches, France 8-18 June 1999. 



1 



CONTENTS 



SYNOPSIS 3 

1. Measuring Higgs boson couplings at the LHC 4 

D. Zeppenfeld, R. Kinnunen, A. Nikitenko and E. Richter-W§s. 

2. Higgs boson production at hadron colliders at NLO 20 

C. Balazs, A. Djouadi, V. Ilyin and M. Spira. 

3. Signatures of Heavy Chcirged Higgs Bosons at the LHC 36 

K.A. Assamagan, A. Djouadi, M. Drees. M. Guchait, R. Kinnunen, J.L. Kneur, 

D. J. Miller, S. Moretti, K. Odagiri and D.R Roy. 

4. Light stop effects and Higgs boson searches at the LHC. 54 

G. Belanger, F. Boudjema, A. Djouadi, V. Ilyin, J.L. Kneur, S. Moretti, 

E. Richter-Wg-s and K. Sridhar. 

5. Double Higgs production at TeV Colliders in the MSSM 67 
R. Lafaye, D.J. Miller, M. Miihlleitner and S. Moretti. 

6. Programs and Tools for Higgs Bosons 88 

E. Boos, A. Djouadi, N. Ghodbane, S. Heinemeyer, V. Ilyin, J. Kalinowski, 
J.L. Kneur and M. Spira. 



2 



SYNOPSIS 



During this Workshop, the Higgs working group has addressed the prospects for searches 
for Higgs particles at future TeV colhders [the Tevatron Runll, the LHC and a future high- 
energy e^e" hnear colhder] in the context of the Standard Model (SM) and its supersym- 
metric extensions such as the Minimal Supersymmetric Standard Model (MSSM). 

In the past two decades, the main focus in Higgs physics at these colliders was on the 
assessment of the discovery of Higgs particles in the simplest experimental detection channels. 
A formidable effort has been devoted to address this key issue, and there is now little doubt 
that a Higgs particle in both the SM and the MSSM cannot escape detection at the LHC or 
at the planed TeV linear e+e~ colhders. 

Once Higgs particles will be found, the next important step and challenge would be to 
make a detailed investigation of their fundamental properties and to establish in all its facets 
the electrowcak symmetry breaking mechanism. To undertake this task, more sophisticated 
analyses are needed since for instance, one has to include the higher-order corrections [which 
are known to be rather large at hadron colliders in particular] to the main detection channels 
to perform precision measurements and to consider more complex Higgs production and 
decay mechanisms [for instance the production of Higgs bosons with other particles, leading 
to multi-body final states] to pin down some of the Higgs properties such as the self-coupling 
or the coupling to heavy states. 

We have addressed these issues at the Les Houches Workshop and initiated a few theo- 
retical/experimental analyses dealing with the measurement of Higgs boson properties and 
higher order corrections and processes. This report summarizes our work. 

The first part of this report deals with the measurements at the LHC of the SM Higgs 
boson couplings to the gauge bosons and heavy quarks. In part 2, the production of the SM 
and MSSM neutral Higgs bosons at hadron colhders, including the next-to-leading order 
QCD radiative corrections, is discussed. In part 3, the signatures of heavy charged Higgs 
particles in the MSSM are analyzed at the LHC. In part 4, the effects of light top squarks 
with large mixing on the search of the lightest MSSM Higgs boson is analyzed at the LHC. 
In part 5, the double Higgs production is studied at hadron and e'^e" colliders in order to 
measure the trihnear Higgs couplings and to reconstruct the scalar potential of the MSSM. 
Finally, part 6 summarizes the work performed on the programs and tools which allow the 
determination of the Higgs boson decay modes and production cross sections at various 
colliders. 

Acknowledgements: 

We thank the organizers of this Workshop, and in particular "le Grand Ordonateur" Patrick 
Aurenche, for the warm, friendly and very stimulating atmosphere of the meeting. We thank 
also our colleagues of the QCD and SUSY working groups for the nice and stimulating, strong 
and super, interactions that we had. Thanks also go to the "personnel" of the Les Houches 
school for allowing us to do physics late at night and for providing us with a hospitable 
environment for many hot or relaxed discussions. 
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Measuring Higgs boson couplings at the LHC 



D. Zeppenfeld, R. Kinnunen, A. Nikitenko and E. Richter-Was 



Abstract 



For an intermediate mass Higgs boson with SM-like couplings the LHC allows 
observation of a variety of decay channels in production by gluon fusion and weak boson 
fusion. Cross section ratios provide measurements of various ratios of Higgs couplings, 
with accuracies of order 15% for 100 fb~^ of data in each of the two LHC experiments. 
For Higgs masses above 120 GeV, minimal assumptions on the Higgs sector allow for 
an indirect measurement of the total Higgs boson width with an accuracy of 10 to 20%, 
and of the H — > WW partial width with an accuracy of about 10%. 

1 Introduction 

Investigation of the symmetry breaking mechanism of the electroweak SU{2) x U{1) gauge 
symmetry will be one of the prime tasks of the LHC. Correspondingly, major efforts have 
been concentrated on devising methods for Higgs boson discovery, for the entire mass range 
allowed within the Standard Model (SM) (100 GeV < < 1 TeV, after LEP2), and for 
Higgs boson search in extensions of the SM, like its minimal supersymmetric extension the 
MSSM |]1], 0. While observation of one or more Higgs scalar(s) at the LHC appears assured, 
discovery will be followed by a more demanding task: the systematic investigation of Higgs 
boson properties. Beyond observation of the various CP even and CP odd scalars which 
nature may have in store for us, this means the determination of the couplings of the Higgs 
boson to the known fermions and gauge bosons, i.e. the measurement of Htt, Hbb, Htt and 
HWW, HZZ, Hj'j couplings, to the extent possible. 

Clearly this task very much depends on the expected Higgs boson mass. For rriH > 
200 GeV and within the SM, only the H ZZ and H WW channels are expected to be 
observable, and the two gauge boson modes are related by SU(2). Above tjih 250 GeV, 
where detector effects will no longer dominate the mass resolution of the H — >■ ZZ — >■ Ai 
resonance, additional information is expected from a direct measurement of the total Higgs 
boson width, Th- A much richer spectrum of decay modes is predicted for the intermediate 
mass range, i.e. if a SM-like Higgs boson has a mass between the reach of LEP2 ( < 110 GeV) 
and the Z-pair threshold. The main reasons for focusing on this range are present indications 
from electroweak precision data, which favor mn < 250 GeV |^, as well as expectations 
within the MSSM, which predicts the lightest Higgs boson to have a mass rrih < 130 GeV @]. 

Until recently, the prospects of detailed and model independent coupling measurements 
at the LHC were considered somewhat remote [§, because few promising search channels 
were known to be accessible, for any given Higgs boson mass. Taking ATLAS search scenarios 
as an example, these were [jl| 



gg ^ H ^ , 
gg^H^ZZ*^M, 



for mj^ < 150 GeV , 
for mH> ISO GeV , 



(1) 
(2) 



4 



and 

gg^H ^ WW* ivlv , for mn > 150 GeV , (3) 

with the possibihty of obtaining some additional information from processes hke WH and/or 
ttH associated production with subsequent H bb and ^ 77 decay for Higgs boson 
masses near 100 GeV. Throughout this contribution, ^^gg H" stands for inclusive Higgs 
production, which is dominated by the gluon fusion process for a SM-like Higgs boson. 

This relatively pessimistic outlook is changing considerably now, due to the demonstra- 
tion that weak boson fusion is a promising Higgs production channel also in the interme- 
diate mass range. Previously, this channel had only been explored for Higgs masses above 
300 GeV. Specifically, it was recently shown in parton level analyses that the weak boson 
fusion channels, with subsequent Higgs decay into photon pairs 0, 0] , 

qq qqH, — 77 , for rriH < 150 GeV , (4) 

into r+r" pairs [0, |, |], 

qq qqH, H tt , for uih < 140 GeV , (5) 

or into W pairs [0, 



qq 



qqH, H WW^*^ e^fi^^r , for > 120 GeV , (6) 



can be isolated at the LHC. Preliminary analyses, which try to extend these parton level 
results to full detector simulations, look promising [^. The weak boson fusion channels 
utilize the significant background reductions which are expected from double forward jet 
tagging 1^ [1^, [l^ and central jet vetoing techniques [|^, |16|, and promise low background 



environments in which Higgs decays can be studied in detail. The parton level results predict 
highly significant signals with (substantially) less than 100 fb^^. 

The prospect of observing several Higgs production and decay channels, over the entire 
intermediate mass range, suggests a reanalysis of coupling determinations at the LHC |^. 
This contribution attempts a first such analysis, for the case where the branching fractions of 
an intermediate mass Higgs resonance are fairly similar to the SM case, i.e. we analyze a SM- 
like Higgs boson only. We make use of the previously published analyses for the inclusive 
Higgs production channels |] and of the weak boson fusion channels P, [fj |8|, P, p!0[|. 
The former were obtained by the experimental collaborations and include detailed detector 
simulations. The latter are based on parton level results, which employ full QCD tree level 
matrix elements for all signal and background processes. We will not discuss here differences 
in the performance expected for the ATLAS and CMS detectors nor details in the theoretical 
assumptions which lead to different estimates for expected signal and background rates. 
The reader is referred to the original publications from which numbers are extracted. In 
Section |^ we summarize expectations for the various channels, including expected accuracies 
for cross section measurement of the various signals for an integrated luminosity of 100 fb~^. 
Implications for the determination of coupling ratios and the measurement of Higgs boson 
(partial) decay widths are then obtained in Section A final summary is given in Section ^ 
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2 Survey of intermediate mass Higgs channels 



The various Higgs channels hsted in Eqs. (|l|-P) and their observabihty at the LHC have all 
been discussed in the literature. Where available, we give values as presently quoted by the 
experimental collaborations. In order to compare the accuracy with which the cross sections 
of different Higgs production and decay channels can be measured, we need to unify these 
results. For example, ii'-factors of unity are assumed throughout. Our goal in this section 
is to obtain reasonable estimates for the relative errors, Ao"h/o"h, which are expected after 
collecting 100 fb~^ in each the ATLAS and the CMS detector, i.e. we estimate results after 
a total of 200 fb~^ of data have been collected at the LHC. Presumably these data will be 
taken with a mix of both low and high luminosity running. 

Table 1: Number of expected events for the inclusive SM if 77 signal and expected 
backgrounds, assuming an integrated luminosity of 100 fb~^ and high luminosity perfor- 
mance. Numbers correspond to optimal 77 invariant mass windows for CMS and ATLAS. 
The expected relative statistical errors on the signal cross section are given for the individual 
experiments and are combined in the last line. 





rriH 


100 


110 


120 


130 


140 


150 


CMS 




18 




Ns 


865 


1038 


1046 


986 


816 


557 










Nb 


29120 


22260 


16690 


12410 


9430 


7790 










AaH/(TH 


20.0% 


14.7% 


12.7% 


11.7% 


12.4% 


16.4% 


ATLAS 1 


i 




Ns 


1045 


1207 


1283 


1186 


973 


652 










Nb 


56450 


47300 


39400 


33700 


28250 


23350 












22.9% 


18.2% 


15.7% 


15.7% 


17.6% 


23.8% 


Combined 


Acfh/cth 


15.1% 


11.4% 


9.9% 


9.4% 


10.1% 


13.5% 



We find that the measurements are largely dominated by statistical errors. For all chan- 
nels, event rates with 200 fb~^ of data will be large enough to use the Gaussian approxima- 
tion for statistical errors. The experiments measure the signal cross section by separately 
determining the combined signal -|- background rate, Ns+b, and the expected number of 
background events, {Nb}- The signal cross section is then given by 

^ Ns+B - {Nb) ^ Ns 

eJCdt eJCdt ' ^ ' 

where e denotes efficiency factors. Thus the statistical error is given by 

^ y/N^ ^ VNs + Nb , . 

an Ns Ns ' ^ ^ 

where in the last step we have dropped the distinction between the expected and the ac- 
tual number of background events. Systematic errors on the background rate are added in 
quadrature to the background statistical error, \/Nb, where appropriate. 
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Well below the H WW threshold, the search for — 77 events is arguably the 
cleanest channel for Higgs discovery. LHC detectors have been designed for excellent two- 
photon invariant mass resolution, with this Higgs signal in mind. We directly take the 
expected signal and background rates for the inclusive if — > 77 search from the detailed 
studies of the CMS and ATLAS collaborations ||T^, which were performed for an 

integrated luminosity of 100 fb~^ in each detector. Expectations are summarized in Table ^ 
Rates correspond to not including a i^'-factor for the expected signal and background cross 
sections in CMS and ATLAS. Cross sections have been determined with the set MRS (Rl) 
of parton distribution functions (pdf's) for CMS, while ATLAS numbers are based on the 
set CTEQ2L of pdf's. 

The inclusive if — 77 signal will be observed as a narrow 77 invariant mass peak on 
top of a smooth background distribution. This means that the background can be directly 
measured from the very high statistics background distribution in the sidebands. We expect 
any systematic errors on the extraction of the signal event rate to be negligible compared to 
the statistical errors which are given in the last row of Table |l[ With 100 fb"^ of data per 
experiment (y{gg —>■ H) ■ B{H —>■ 77) can be determined with a relative error of 10 to 15% 
for Higgs masses between 100 and 150 GeV. Here we do not include additional systematic 
errors, e.g. from the luminosity uncertainty or from higher order QCD corrections, because 
we will mainly consider cross section ratios in the final analysis in the next Section. These 
systematic errors largely cancel in the cross section ratios. Systematic errors common to 
several channels will be considered later, where appropriate. 

A Higgs search channel with a much better signal to background ratio, at the price of 
lower statistics, however, is available via the inclusive search for H ZZ* — >■ Ai events. 
Expected event numbers for 100 fb"^ in both ATLAS and CMS |l| are listed in Table ^ 
These numbers were derived using CTEQ2L pdf's and are corrected to contain no QCD 
K-factor. For those Higgs masses where no ATLAS or CMS prediction is available, we 
interpolate/extrapolate the results for the nearest Higgs mass, taking the expected H — > ZZ* 
branching ratios into account for the signal. Similar to the case of ii —> 77 events, the signal 
is seen as a narrow peak in the four-lepton invariant mass distribution, i.e. the background 
can be extracted directly from the signal sidebands. The combined relative error on the 
measurement of (j{gg H) ■ B{H ZZ*) is listed in the last line of Table ^. For Higgs 
masses in the 130-150 GeV range, and above Z-pair threshold, a 10% statistical error on 
the cross section measurement is possible. In the intermediate range, where H — > WW 
dominates, and for lower Higgs masses, where the Higgs is expected to dominantly decay 
into bb, the error increases substantially. 

Above niH ^ 135 GeV, H — * becomes the dominant SM Higgs decay channel. 

The resulting inclusive WW — > signal is visible above backgrounds, after exploiting 

the characteristic lepton angular correlations for spin zero decay into W pairs near thresh- 
old [^. The inclusive channel, which is dominated hj gg —y H —y WW, has been analyzed 
by ATLAS for uih > 150 GeV and for integrated luminosities of 30 and 100 fb~^ ||l[ and 
by CMS for uih > 120 GeV and 30 fb^^ PO]. The expected event numbers for 30 fb~^ are 



listed in Table y. The numbers are derived without QCD K-factors and use CTEQ2L for 
ATLAS and MRS (A) pdf's for CMS results. 
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Table 2: Number of expected events for the inclusive SM H ZZ* i^i^i^i^ signal and 
expected backgrounds, assuming an integrated luminosity of 100 fb~^ and high luminosity 
performance. Numbers correspond to optimal four-lepton invariant mass windows for CMS 
and ATLAS and to the combined total. Rates in parentheses correspond to numbers in- 
terpolated, according to if — > ZZ* branching ratios for the signal. The expected relative 
statistical errors on the signal cross section are given for each experiment and are combined 
in the last line. 





120 


130 


140 


150 


160 


170 


180 


CMS 19 


Ns 


19.2 


55.3 


(99) 


131.4 


(48) 


29.4 


(76.5) 




Nb 


12.9 


17.1 


(20) 


22.5 


(26) 


27.5 


(27) 






29.5% 


15.4% 


11.0% 


9.4% 


17.9% 


25.7% 


13.3% 


ATLAS i 


Ns 


10.3 


28.7 


(51) 


67.6 


(31) 


19.1 


49.7 




Nb 


4.44 


7.76 


(8) 


8.92 


(8) 


8.87 


8.81 






37.3% 


21.0% 


15.1% 


12.9% 


20.1% 


27.7% 


15.4% 


Combined 




23.1% 


12.4% 


8.9% 


7.6% 


13.4% 


18.8% 


10.1% 



Unlike the two previous modes, the two missing neutrinos in the B. — >■ WW events do 
not allow for a reconstruction of the narrow Higgs mass peak. Since the Higgs signal is only 
seen as a broad enhancement of the expected background rate in lepton-neutrino transverse 
mass distributions, with similar shapes of signal and background after application of all cuts, 
a precise determination of the background rate from the data is not possible. Rather one 
has to rely on background measurements in phase space regions where the signal is weak, 
and extrapolation to the search region using NLO QCD predictions. The precise error on 
this extrapolation is unknown at present, the assumption of a 5% systematic background 
uncertainty appears optimistic but attainable. It turns out that with 30 fb~^ already, the 
systematic error starts to dominate, because the background exceeds the signal rate by fac- 
tors of up to 5, depending on the Higgs mass. Running at high luminosity makes matters 
worse, because the less efficient reduction of it backgrounds, due to less stringent 6-jet veto 
criteria, increases the background rate further. Because of this problem we only present 
results for 30 fb~^ of low luminosity running in Table |^. Since neither of the LHC collabo- 
rations has presented predictions for the entire Higgs mass range, we take CMS simulations 
below 150 GeV and ATLAS results at 190 GeV, but divide the resultant statistical errors by a 
factor a/2, to take account of the presence of two experiments. Between 150 and 180 GeV we 
combine both experiments, assuming 100% correlation in the systematic 5% normalization 
error of the background. 

The previous analyses are geared towards measurement of the inclusive Higgs production 
cross section, which is is dominated by the gluon fusion process. 15 to 20% of the signal 
sample, however, is expected to arise from weak boson fusion, qq qqH or corresponding 
antiquark initiated processes. The weak boson fusion component can be isolated by making 
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Table 3: Number of expected events for the inclusive SM H WW* V^vi~v signal and 
expected backgrounds, assuming an integrated luminosity of 30 fb~^. Numbers correspond to 
optimized cuts, varying with the mass of the Higgs boson being searched for. The expected 
relative errors on the signal cross section are given for each experiment, separating the 
statistical error, the effect of a systematic 5% error of the background level, and the two 
added in quadrature. The combined error for the two experiments assumes 100% correlation 
of the systematic errors on the background determination. 



rau 


120 


130 


140 


150 


160 


170 


180 


190 


CMS 


Ns 


44 


106 


279 


330 


468 


371 


545 




m 


Nb 


272 


440 


825 


732 


360 


360 


1653 






A(j/f/(T_f/(stat.) 


40.4% 


22.0% 


11.9% 


9.9% 


6.1% 


7.3% 


8.6% 






Ao-///o-//(syst.) 


30.9% 


20.8% 


14.8% 


11.1% 


3.8% 


4.9% 


15.2% 






A(7i^ /(Tj^ (comb . ) 


50.9% 


30.3% 


19.0% 


14.9% 


7.3% 


8.8% 


17.4% 


20.6% 


ATLAS 


Ns 








240 


400 


337 


276 


124 


i 


Nb 








844 


656 


484 


529 


301 


Aan/f^H (stat.) 








13.7% 


8.1% 


8.5% 


10.3% 


16.6% 




Aan/cTH (syst.) 








17.6% 


8.2% 


7.2% 


9.6% 


12.1% 




Aan/o'H (comb.) 


50.9% 


30.3% 


19.0% 


22.3% 


11.5% 


11.1% 


14.1% 


20.6% 


Com 


Aan/cTH (comb.) 


42.1% 


26.0% 


17.0% 


14.8% 


7.0% 


8.0% 


13.6% 


16.9% 



use of the two forward tagging jets which are present in these events and by vetoing additional 
central jets, which are unlikely to arise in the color singlet signal process A more detailed 
discussion of these processes can be found in Ref. [0 from which most of the following 
numbers are taken. 

The qq — * qqH, H ^ 'jj process was first analyzed in Ref. P], where cross sections for 
signal and background were obtained with full QCD tree level matrix elements. The parton 
level Monte Carlo determines all geometrical acceptance corrections. Additional detector 
effects were included by smearing parton and photon 4-momenta with expected detector 
resolutions and by assuming trigger, identification and reconstruction efficiencies of 0.86 for 
each of the two tagging jets and 0.8 for each photon. Resulting cross sections were presented 
in Ref. for a fixed 77 invariant mass window of total width Am^^ = 2 GeV. We correct 
these numbers for mn dependent mass resolutions in the experiments. We take lAa mass 
windows, as given in Ref. IQ for high luminosity running, which are expected to contain 79% 
of the signal events for ATLAS. The 2 GeV window for tuh = 100 GeV at CMS 0, |l8l is 



assumed to scale up like the ATLAS resolution and assumed to contain 70% of the Higgs 
signal. The expected total signal and background rates for 100 fb~^ and resulting relative 
errors for the extraction of the signal cross section are given in Table ^. Statistical errors only 
are considered for the background subtraction, since the background level can be measured 
independently by considering the sidebands to the Higgs boson peak. 

The next weak boson fusion channel to be considered is qq — >■ qqH, H tt. Again, 
this channel has been analyzed at the parton level, including some estimates of detector 
effects, as discussed for the if — >• 77 case. Here, a lepton identification efficiency of 0.95 
is assumed for each lepton ^ = e, /i. Two r-decay modes have been considered so far: 
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Table 4: Number of expected 77jj events from the qq qqH, — 77 weak boson 
fusion signal and expected backgrounds, assuming an integrated luminosity of 100 fb~^. 
Numbers correspond to optimal 77 invariant mass windows for CMS and ATLAS and to 
the combined total, as projected from the parton level analysis of Refs. ^. The expected 
relative statistical errors on the signal cross section are given for each experiment and are 
combined in the last line. 





rriH 


100 


110 


120 


130 


140 


150 


projected CMS 


Ns 


37 


48 


56 


56 


48 


33 


performance 


Nb 


33 


32 


31 


30 


28 


25 






22.6% 


18.6% 


16.7% 


16.6% 


18.2% 


23.1% 


projected ATLAS 


Ns 


42 


54 


63 


63 


54 


37 


performance 


Nb 


61 


60 


56 


54 


51 


46 






24.2% 


19.8% 


17.3% 


17.2% 


19.0% 


24.6% 


combined 




16.5% 


13.6% 


12.0% 


11.9% 


13.1% 


16.8% 



Table 5: Number of expected signal and background events for the qq — > qqH TTjj 
channel, for 100 fb~^ and two detectors. Cross sections are added for rr (!^h^;^x and 
rr e^fi^^T events as given in Refs. |]^, 0. The last line gives the expected statistical 
relative error on the qq — > qqH, H ^ tt cross section. 



rriH 


100 


110 


120 


130 


140 


150 


Ns 


211 


197 


169 


128 


79 


38 


Nb 


305 


127 


51 


32 


27 


24 


^cth/c^h 


10.8% 


9.1% 


8.8% 


9.9% 


13.0% 


20.7% 



if — >• rr ^ i'^h^'^T lH and H ^ tt ^ e^fi^;^T |H- These analyses were performed for low 
luminosity running. Some deterioration at high luminosity is expected, as in the analogous 
H/A TT channel in the MSSM search [jl|. At high luminosity, pile- up effects degrade 
the '^T resolution significantly, which results in a worse rr invariant mass resolution. At 
a less significant level, a higher px threshold for the minijet veto technique will increase 
the QCD and ti backgrounds. The r-identification efficiency is similar at high and low 
luminosity. We expect that the reduced performance at high luminosity can be compensated 
for by considering the additional channels H ^ tt ^ e+e^l^r? A*^/^ I^t- 2'+jets and 
ZZ+jets backgrounds (with ZZ l^l'vv) are strongly suppressed by rejecting same flavor 
lepton pairs which are compatible with Z decays (m« = niz ± 6 GeV). Drell-Yan plus jets 
backgrounds are further reduced by requiring significant Since these analyses have not 
yet been performed, we use the predicted cross sections for only those two channels which 
have already been discussed in the hterature and scale event rates to a combined 200 fb~^ 
of data. Results are given in Table ^ 
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Table 6: Number of events expected for qq — > qqH, H WW^*'^ /i^e^l^T in 200 fb~^ 
of data, and corresponding backgrounds The expected relative statistical error on the 
signal cross section is given in the last line. 



The previous two weak boson channels allow reconstruction of the Higgs resonance as 
an invariant mass peak. This is not the case for H —>■ WW — ^ as discussed previ- 

ously for the inclusive search. The weak boson fusion channel can be isolated separately by 
employing forward jet tagging and color singlet exchange isolation techniques in addition to 
tools like charged lepton angular correlations which are used for the inclusive channel. The 
corresponding parton level analysis for qq — > qqH, H — »• — > jjL^e^^T has been per- 
formed in Ref. 0] and we here scale the results to a total integrated luminosity of 200 fb~^, 
which takes into account the availability of two detectors. As for the tau case, the analysis 
was done for low luminosity running conditions and somewhat higher backgrounds are ex- 
pected at high luminosity. On the other hand the — > and e^e~j)T 
modes should roughly double the available statistics since very few signal events have lepton 
pair invariant masses compatible with Z ^ it decays. Therefore our estimates are actually 
conservative. Note that the expected background for this weak boson fusion process is much 
smaller than for the corresponding inclusive measurement. As a result modest systematic 
uncertainties will not degrade the accuracy with which o"(gg — > qqH^ ■ B(H WW^*^) can 
be measured. A 10% systematic error on the background, double the error assumed in the 
inclusive case, would degrade the statistical accuracy by, typically, a factor 1.2 or less. As 
a result, we expect that a very precise measurement of a{qq — > qqH) ■ B{H — > H^iy*^*)) can 
be performed at the LHC, with a statistical accuracy of order 5% or even better in the mass 
range rriH > 140 GeV. Even for rriH as low as 120 GeV a 12% measurement is expected. 

3 Measurement of Higgs properties 

One would like to translate the cross section measurements of the various Higgs production 
and decay channels into measurements of Higgs boson properties, in particular into measure- 
ments of the various Higgs boson couplings to gauge fields and fermions. This translation 
requires knowledge of NLO QCD corrections to production cross sections, information on 
the total Higgs decay width and a combination of the measurements discussed previously. 
The task here is to find a strategy for combining the anticipated LHC data without undue 
loss of precision due to theoretical uncertainties and systematic errors. 

For our further discussion it is convenient to rewrite all Higgs boson couplings in terms 
of partial widths of various Higgs boson decay channels. The Higgs-fermion couplings gnff, 
for example, which in the SM are given by the fermion masses, gnff = ^fi^H)/^, can be 



rriH 



120 130 140 150 160 170 180 190 



Ns 
Nb 



136 332 592 908 1460 1436 1172 832 
136 160 188 216 240 288 300 324 
12.1% 6.7% 4.7% 3.7% 2.8% 2.9% 3.3% 4.1% 
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traded for the H ^ f f partial widths, 

Here c/ is the color factor (1 for leptons, 3 for quarks). Similarly the square of the HWW 
coupling {gnww = gf^w in the SM) or the HZZ coupling is proportional to the partial 
widths Tw = r(if — > lyiy*) or Tz = T{H — > ZZ*) [^. Analogously we trade the squares 



of the effective H^yy and Hgg couplings for = T{H ^ 77) and Tg = T{H ^ gg). Note 
that the Hgg coupling is essentially proportional to gnu, the Higgs boson coupling to the 
top quark. 

The Higgs production cross sections are governed by the same squares of couplings. This 



allows to write e.g. the gg ^ H production cross section as p2 | 

(^{gg ^ H) = V{H ^ gg)—^T —g{x,m]j)g{-,m]j), (10) 

where r = m'j^/s. Similarly the qq qqH cross sections via WW and ZZ fusion are 
proportional to r[H —>■ WW*) and r{H ZZ*), respectively. In the narrow width ap- 
proximation, which is appropriate for the intermediate Higgs mass range considered here, 
these production cross sections need to be multiplied by the branching fractions for final 
state j, B{H ^ j) = Tj/T, where T denotes the total Higgs width. This means that the 
various cross section measurements discussed in the previous Section provide measurements 
of various combinations TiTj/T. 

The production cross sections are subject to QCD corrections, which introduces theo- 



retical uncertainties. While the /T-factor for the gluon fusion process is large [^], which 
suggests a sizable theoretical uncertainty on the production cross section, the NLO correc- 
tions to the weak boson fusion cross section are essentially identical to the ones encountered 
in deep inelastic scattering and are quite small [^. Thus we can assign a small theoretical 
uncertainty to the latter, of order 5%, while we shall use a larger theoretical error for the 



gluon fusion process, of order 20% ||2^ . The problem for weak boson fusion is that it consists 
of a mixture of ZZ H and WW H events, and we cannot distinguish between the 
two experimentally. In a large class of models the ratio of HWW and HZZ couplings is 
identical to the one in the SM, however, and this includes the MSSM. We therefore make 
the following W, Z-universality assumption: 

• The H ZZ* and H WW* partial widths are related by SU(2) as in the SM, i.e. 
their ratio, z, is given by the SM value, 

Tz = z Tw = Tw ■ (11) 



Note that this assumption can be tested, at the 15-20% level for uih > 130 GeV, by forming 
the ratio Ba{gg ^ H ^ ZZ*)/Ba{gg ^ H ^ WW*), in which QCD uncertainties cancel 
(see Table ^. 
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With W, Z-universality, the three weak boson fusion cross sections give us direct mea- 
surements of three combinations of (partial) widths, 

r r 

= — ^ from qq — > qqH, — > 77 , (12) 

Xr = — from qq qqH, H ^ tt , (13) 

Xw = ^ from qq qqH, H WW^*^ , (14) 

with common theoretical systematic errors of 5%. In addition the three gluon fusion channels 
provide measurements of 

"p "p 

(15) 
(16) 
(17) 





r r 
r 


from 


99 - 


^ H - 




Yz = 


r 


from 


99 - 


^ H - 




Yw = 


r 


from 


99 - 


H - 





with common theoretical systematic errors of 20%. 

The first precision test of the Higgs sector is provided by taking ratios of the Xj's and 
ratios of the YiS. In these ratios the QCD uncertainties, and all other uncertainties related 
to the initial state, like luminosity and pdf errors, cancel. Beyond testing W, Z-universality, 
these ratios provide useful information for Higgs masses between 100 and 150 GeV and 120 to 
150 GeV, respectively, where more than one channel can be observed in the weak boson fusion 
and gluon fusion groups. Typical errors on these cross section ratios are expected to be in 
the 15 to 20% range (see Table 0). Accepting an additional systematic error of about 20%, a 
measurement of the ratio Tg/T\Y, which determines the Htt to HWW coupling ratio, can be 
performed, by measuring the cross section ratios Ba{gg — if — >• 77)/cr(gg qqH)B{H — > 
77) and Ba{gg ^ H ^ WW*)/a{qq qqH)B{H WW*). Expected accuracies are 
listed in Table 0. In these estimates the systematics coming from understanding detector 
acceptance is not included. 

Beyond the measurement of coupling ratios, minimal additional assumptions allow an 
indirect measurement of the total Higgs width. First of all, the r partial width, properly 
normalized, is measurable with an accuracy of order 10%. The r is a third generation 
fermion with isospin — |, just like the 6-quark. In all extensions of the SM with a common 
source of lepton and quark masses, even if generational symmetry is broken, the ratio of h 
to r Yukawa couplings is given by the fermion mass ratio. We thus assume, in addition to 
W, Z-universality, that 

• The ratio of 6 to r couplings of the Higgs is given by their mass ratio, i.e. 

y=Tr = ^cqcd— — = ^cqcd 5 — , (18) 

where cqcd is the known QCD and phase space correction factor. 
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Table 7: Summary of the accuracy with which various ratios of partial widths can be deter- 
mined with 200 fb~^ of data. The first two columns give the ratio considered and indicate 
the method by which it is measured. Y^/IV, for example, indicates a measurement of 
aB[H ZZ*)/aB[H — > lyVT*) in gluon fusion, while Xi ratios correspond to weak boson 
fusion (see text for details). The statistical combination of several channels for a given width 
ratio is indicated by ©. 5% and 20% theoretical uncertainties for weak boson and gluon fu- 
sion cross sections affect the mixed gluon/weak boson fusion ratios only, which are needed 
for a measurement of Tg/Tiy. The effect of this systematic error is indicated in the last line. 



rriH 




100 110 120 130 140 150 160 170 180 


z = Tz/Tw 


Yz/Yw 

Yz ^7 

Xz_ Yz ^7 

Y\Y Y'Y 


48% 29% 19% 17% 15% 20% 17% 

30% 21% 19% 23% 

29% 19% 15% 14% 15% 20% 17% 




^7 

Y^ Xw 


16% 12% 11% 13% 
15% 12% 14% 21% 
20% 16% 15% 16% 18% 27% 




1^ © ^ © 21% 


22% 18% 15% 13% 12% 13% 8% 9% 14% 
30% 27% 25% 24% 24% 24% 22% 22% 25% 



• The total Higgs width is dominated by decays to 66, rr, WW, ZZ, gg and 77, i.e. the 
branching ratio for unexpected channels is small: 

e = 1 - B{H bb) + B{H rr) + B{H WW'-*^) + 

B{H ^ ZZ^*'>) + B{H ^gg) + B{H ^-f-f)^ <1. (19) 

Note that, in the Higgs mass range of interest, these two assumptions are satisfied for both 
CP even Higgs bosons in most of the MSSM parameter space. The first assumption holds in 
the MSSM at tree level, but can be violated by large squark loop contributions, in particular 
for small and large tan /3 ^ . The second assumption might be violated, for example, 
if the H ^ cc partial width is exceptionally large. However, a large up-type Yukawa coupling 
would be noticeable in the Tg/Tw coupling ratio, which measures the Htt coupling. 
With these assumptions consider the observable 

Tw = Xr{l + y)+Xwil + z)+X^ + Xg 

= (Tr + Tt + Tw + Tz + T^ + Tg^^ = {l-e)Tw , (20) 

where Xg = TgTw/^ is determined by combining Yw and the product Y^Xw/X-y. Tw 
provides a lower bound on T{H — > lyiy^*)) = T]y. Provided e is small (e < 0.1 suffices for 
practical purposes), the determination of Tw provides a direct measurement of the H 
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Vrvr^*^ partial width. Once Tw has been determined, the total width of the Higgs boson is 
given by 

^ ^Xr{l + y)+Xw{l + z)+X^ + Xg) \, ^ (21) 



r 



J- w 



x^ x^V"'' ■ ■ ■ ^ V (1-6)2 ■ 

For a SM-like Higgs boson the Higgs width is dominated by the H ^bb and H WW^*^ 
channels. Thus, the error on Tw is dominated by the uncertainties of the X^y and X^. 
measurements and by the theoretical uncertainty on the 6-quark mass, which enters the 
determination of y quadratically. According to the Particle Data Group, the present uncer- 
tainty on the b quark mass is about ±3.5% ||2^. Assuming a luminosity error of ±5% in 



addition to the theoretical uncertainty of the weak boson fusion cross section of ±5%, the 
statistical errors of the qq — > qqH, H ^ tt and qq qqH, H WW cross sections of 
Tables ^ and |^ lead to an expected accuracy of the Tw determination of order 10%. More 
precise estimates, as a function of the Higgs boson mass, are shown in Fig. |I|. 
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Figure 1: Expected accuracy with which the Higgs boson width can be measured at the 
LHC, with 100 fb~^ of data in each experiment. Results are shown for the extraction of the 
the H — i> WW partial width, Tw, and and the total Higgs boson width, F. e is the sum of 
the residual (small) branching ratios of unobserved channels, mainly H ^ cc (see text). 



The extraction of the total Higgs width, via Eq. (|21|), requires a measurement of the qq - 
qqH, H — > cross section, which is expected to be available for rriH > 115 GeV [|10 

Consequently, errors are large for Higgs masses close to this lower limit (we expect a relative 
error of ~ 20% for mu = 120 GeV and e < 0.05). But for Higgs boson masses around the 
WW threshold, F(l — e)^ can be determined with an error of about 10%. Results are shown 
in Fig. and look highly promising. 
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4 Summary 



In the last section we have found that various ratios of Higgs partial widths can be measured 
with accuracies of order 10 to 20%, with an integrated luminosity of 100 fb~^ per experiment. 
This translates into 5 to 10% measurements of various ratios of couphng constants. The ratio 
Ft/Fw measures the coupling of down-type fermions relative to the Higgs couplings to gauge 
bosons. To the extent that the H'-fj triangle diagrams are dominated by the W loop, the 
width ratio Tr/T-y measures the same relationship. The fermion triangles leading to an 
effective Hgg coupling are expected to be dominated by the top-quark, thus, Tg/Tw probes 
the couphng of up-type fermions relative to the HWW coupling. Finally, for Higgs boson 
masses above ^ 120 GeV, the absolute normalization of the HWW coupling is accessible 
via the extraction of the H WW^*'' partial width in weak boson fusion. 

Note that these measurements test the crucial aspects of the Higgs sector. The HWW 
coupling, being linear in the Higgs field, identifies the observed Higgs boson as the scalar 
responsible for the spontaneous breaking of SU{2) x U{1): a scalar without a vacuum ex- 
pectation value couples to gauge bosons only via HHWW or HHW vertices at tree level, 
i.e. the interaction is quadratic in scalar fields. The absolute value of the HWW coupling, 
as compared to the SM expectation, reveals whether H may be the only mediator of spon- 
taneous symmetry breaking or whether additional Higgs bosons await discovery. Within the 
framework of the MSSM this is a measurement of | sin(/3 — q;)|, at the ±0.05 level. The 
measurement of the ratios of gmt/ 9hww ^-nd Qhtt/ Qeww then probes the mass generation 
of both up and down type fermions. 

The results presented here constitute a first look only at the issue of coupling extractions 
for the Higgs. This is the case for the weak boson fusion processes in particular, which prove 
to be extremely valuable if not essential. Our analysis is mostly an estimate of statistical 
errors, with some rough estimates of the systematic errors which are to be expected for the 
various measurements of (partial) widths and their ratios. A number of issues need to be 
addressed in further studies, in particular with regard to the weak boson fusion channels. 

(a) The weak boson fusion channels and their backgrounds have only been studied at the 
parton level, to date. Full detector level simulations, and optimization of strategies 
with more complete detector information is crucial for further progress. 

(b) A central jet veto has been suggested as a powerful tool to suppress QCD backgrounds 

to the color singlet exchange processes which we call weak boson fusion. The feasibility 
of this tool and its reach need to be investigated in full detector studies, at both low 
and high luminosity. 

(c) In the weak boson fusion studies of i7 — > WW and H ^ tt decays, double leptonic 
e~^e~^T and fi~^fi~^T signatures have not yet been considered. Their inclusion promises 
to almost double the statistics available for the Higgs coupling measurements, at the 
price of additional ZZ+jets and Drell-Yan plus jets backgrounds which are expected 
to be manageable. 
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(d) Other channels, hke WH or tiH associated production with subsequent decay H ^bb 
OT H ^ 77, provide additional information on Higgs coupling ratios, which complement 
our analysis at small Higgs mass values, mn ^ 120 GeV 0, |^. These channels need 
to be included in the analysis. 

(e) Much additional work is needed on more reliable background determinations. For the 
H WW^*^ — i> ~j)T channel in particular, where no narrow Higgs resonance peak 
can be reconstructed, a precise background estimate is crucial for the measurement 
of Higgs couplings. Needed improvements include NLO QCD corrections, single top 
quark production backgrounds, the combination of shower Monte Carlo programs with 
higher order QCD matrix element calculations and more. 

(f ) Both in the inclusive and WBF analyses any given channel contains a mixture of events 
from gg ^ H and qq — > qqH production processes. The determination of this mixture 
adds another source of systematic uncertainty, which was not included in the present 
study. In ratios of X observables (or of different Yi) these uncertainties largely cancel, 
except for the effects of acceptance variations due to different signal selections. Since 
an admixture from the wrong production channel is expected at the 10 to 20% level 
only, these systematic errors are not expected to be serious. 

(g) We have only analyzed the case of a single neutral, CP even Higgs resonance with 
couplings which are close to the ones predicted in the SM. While this case has many 
applications, e.g. for the large region of the MSSM, more general analyses, in 
particular of the MSSM case, are warranted and highly promising. 

While much additional work is needed, our study clearly shows that the LHC has excellent 
potential to provide detailed and accurate information on Higgs boson interactions. The 
observability of the Higgs boson at the LHC has been clearly established, within the SM and 
extensions like the MSSM. The task now is to sharpen the tools for accurate measurements 
of Higgs boson properties at the LHC. 
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Higgs boson production at hadron colliders at NLO 

C. Balazs, a. Djouadi, V. Ilyin and M. Spira 

Abstract 

We discuss the production of neutral Higgs bosons at the hadron colhders Teva- 
tron and LHC, in the context of the Standard Model and its minimal super symmetric 
extension. The main focus will be on the next-to-leading order QCD radiative correc- 
tions to the main Higgs production mechanisms and on Higgs production in processes 
of higher order in the strong coupling constant. 

1 Introduction 

One of the most important missions of future high-energy colliders will be the search for 
scalar Higgs particles and the exploration of the electroweak symmetry breaking mechanism. 
In the Standard Model (SM), one doublet of complex scalar fields is needed to spontaneously 
break the symmetry, leading to a single neutral Higgs particle [|l|. In the SM, the Higgs 
boson mass is a free parameter and can have a value anywhere between 100 GeV and 1 TeV. 
In contrast, a firm prediction of supersymmetric extensions of the SM is the existence of a 
light scalar Higgs boson In the Minimal Supersymmetric Standard Model (MSSM) the 
Higgs sector contains a quintet of scalar particles [two CP-even h and H, a pseudoscalar A 
and two charged particles] the Higgs boson h of which should be light, with a mass 
Mh < 135 GeV. If this particle is not found at LEP2, it will be produced at the upgraded 
Tevatron (where a large luminosity, / £ ~ 20 fb~^, is expected) [||, |^ or at the LHC 0, ^ |^, 
if the MSSM is indeed realized in Nature. 

Since Higgs boson production at hadron colliders involves strongly interacting particles 
in the initial state, the lowest order cross sections are in general affected by large uncertain- 
ties arising from higher order corrections. If the next-to-leading QCD corrections to these 
processes are included, the total cross sections can be defined properly and in a reUable way 
in most of the cases. In this contribution, we will discuss the next -to-leading order (NLO) 
QCD radiative corrections to the main neutral Higgs production mechanisms as well as 
neutral Higgs boson production in processes of higher order in the strong coupling constant. 

The contribution is organized as follows. In the next section 0, we summarize the main 
processes for the production of the neutral Higgs bosons of the MSSM at hadron colliders 
and discuss the effects of their next-to-leading order QCD corrections; we will then discuss 
the recently evaluated SUSY-QCD corrections to some of these processes. In section 3 
we will concentrate on Higgs boson production in association with heavy quarks which in the 
MSSM might have the largest cross sections due a possible strong enhancement of the Yukawa 
couplings of third generation quarks; we will discuss in particular the next-to-leading order 
QCD corrections to Higgs production in heavy quark fusion. In section 4 §], we will analyze 
the detection of the SM and lightest MSSM [in the decoupling regime] Higgs boson in the 
channel 77+jet at the LHC [where the Higgs boson is produced in the gluon-gluon fusion 
mechanism and decays into two photons]. 
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2 MSSM neutral Higgs production at hadron colliders: 
Next— to— Leading— Order QCD corrections 

2.1 Summary of standard NLO QCD corrections 

At hadron colliders, the production of the neutral Higgs bosons in the MSSM is provided by 
the following processes: 

(a) The gluon-gluon fusion, mediated by heavy quark loops, is the dominant production 
mechanism for neutral Higgs particles, gg ^ with ^ = h, H or A ||10[ . Since the Higgs 
particles in the mass range of interest, Mq, < 135 GeV, dominantly decay into bottom quark 
pairs, this process is rather difficult to exploit at the Tevatron because of the huge QCD 
background 0. In contrast, at the LHC rare decays of the h boson to two photons or decays 
of the H, A bosons to r and fi lepton pairs make this process very useful ||. 

(b) Higgs-strahlung oS W or Z bosons for the CP-even Higgs particles [due to CP- 
invariance, the pseudoscalar A particle does not couple to the massive gauge bosons at tree 
level]: qq V* with ^ = h, H and ^ = |]lll. At the Tevatron, the process 
qq' — >■ hW [with the h boson decaying into bb pairs] develops a cross section of the order of 
a fraction of a picobarn for a SM-like h boson with a mass below ~ 135 GeV, making it the 
most relevant mechanism to study 0. At the LHC, both the bb and 77 decay modes of the 
h boson may be exploited 

(c) If the heavier H, A, bosons are not too massive, the pair production of two Higgs 
particles in the Drell-Yan type process, qq — > $i$2 [P"^H^, might lead to a variety of final 



states [hA, HA, H^h, H^H, H^A, H+H'] with reasonable cross sections [in particular for 
Ma ~ Mh ~ Mh± ^ 250 GeV and small values of tan/3, the ratio of the vacuum expectation 
values of the two Higgs doublets] especially at the LHC. Moreover, neutral and charged Higgs 



boson pairs will be produced in gluon fusion gg — > $i$2 llT3Hl5|. 

(d) The production of CP-even Higgs bosons via vector boson fusion, qq qqV*V* 
qq^ Ijl^. In the case of a SM-like h boson, this process has a sizeable cross section at the 
LHC. While decays of the Higgs boson into heavy quark pairs are problematic to be detected 
in the jetty environment of the LHC, decays into r lepton pairs make this process useful at 
the LHC as discussed recently . 



(e) The production of neutral Higgs bosons via radiation off heavy bottom and top quarks 
[qq, gg —>■ bb^, tt$] might play an important role in SUSY theories [^. In particular, because 
the couplings of the Higgs boson to b quarks can be strongly enhanced for large values of 
tan/3, Higgs production in association with bb pairs can give rise to large production rates. 

It is well known that for processes involving strongly interacting particles, as is the case 
for the ones discussed above, the lowest order cross sections are affected by large uncertain- 
ties arising from higher order corrections. If the next-to-leading QCD corrections to these 
processes are included, the total cross sections can be defined properly and in a rehable way 
in most of the cases. 
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For the standard QCD corrections, the next-to-leading corrections are available for most 
of the Higgs boson production processes^. They are parameterized by the K-factors [defined 
as the ratios of the next-to-leading order cross sections to the lowest order ones]: 

- For Higgs boson production via the gluon fusion processes, the K-factors have been 
calculated a few years ago in the SM |]20| and in the MSSM the [two-loop] QCD 



corrections to the heavy top and to the bottom quark loops [which gives the dominant 
contributions to the cross section for large tan/? values] have been found to be significant 
since they increase the cross sections by up to a factor of two. 

- The K-factors for Higgs production in association with a gauge boson (6) and for Drell- 
Yan-like Higgs pair production (c), can be inferred from the one the Drell-Yan production 
of weak vector bosons and increase the cross section by approximately 30% [^ . 



- The QCD corrections to pair production gg $i$2 are only known in the limit of 
light Higgs bosons compared with the loop-quark mass. This is a good approximation in 
the case of the lightest h boson which, due to phase space, has the largest cross section in 
which the top quark loop is dominant for small values of tan/? or in the decoupling limit. 
The corrections enhance the cross sections by up to a factor of two ||15|| . 

- For Higgs boson production in the weak boson fusion process (rf), the QCD corrections 
can be derived in the structure function approach from deep-inelastic scattering; they turn 
out to be rather small, enhancing the cross section by about 10% [ P^ . 

- Finally, the full QCD corrections to the associated Higgs production with heavy quarks 
(e) are not yet available; they are only known in the limit of light Higgs particles compared 
with the heavy quark mass [Q which is only applicable to tth production; in this limit the 
QCD corrections increase the cross section by about 20-60%. 

2.2 SUSY QCD corrections 

Besides these standard QCD corrections, additional SUSY-QCD corrections must be taken 
into account in SUSY theories; the SUSY partners of quarks and gluons, the squarks and 
gluinos, can be exchanged in the loops and contribute to the next-to-leading order total cross 
sections. In the case of the gluon fusion process, the QCD corrections to the squark loop 
contributions have been calculated in the limit of light Higgs bosons and heavy gluinos; the 
K-factors were found to be of about the same size as the ones for the quark loops ||25|| . 

During this workshop, we studied the SUSY-QCD corrections to the Higgs production 
cross sections for Higgs-strahlung, Drell-Yan like Higgs pair production and weak boson 



fusion processes [^. This analysis completes the theoretical calculation of the NLO pro- 
duction cross sections of these processes in the framework of supersymmetric extensions 
of the Standard Model. These corrections originate from qqV one-loop vertex corrections, 
where squarks of the first two generations and gluinos are exchanged, and the corresponding 
quark self-energy counterterms. Fig. |I[ 



^The small NLO QCD corrections to the important Higgs decays into photons are also available llG 
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X/WW^ V 



Figure 1: Generic diagrams contributing to the SUSY-QCD corrections to the qqV vertex 
\V = 7, Z, W] at next-to-leading order. 



Including these SUSY-particle loop corrections, the lowest order partonic cross section 
for the Drell-Yan type processes will be shifted by 



0"LO o-LO 



TT 



5ReC(s, rriq, nig) 



For degenerate unmixed squarks [as is approximately the case for the first two generation 
squarks], the expression of the factor C is simply given by 



C{s,mg,mg) = 2 / xdx / dylog 



m? + (m? — m?)x 



-sx'^y{l — y) + (wi? — m|)x + m| — ie 

For the fusion processes, the standard QCD corrections have been calculated within 
the structure function approach [Q. Since at lowest order, the proton remnants are color 
singlets, at NLO no color will be exchanged between the first and the second incoming (out- 
going) quark line and hence the QCD corrections only consist of the well-known corrections 
to the structure functions Fj(x,M^) {i = 1,2,3). The final result for the QCD-corrected 
cross section can be obtained from the replacements 



M^) ^ Fiix, M^) + AFi(x, M\ 



1,2,3) 



(3) 



with AFj(x, M^,Q^) the standard QCD corrections The typical renormalization and 
factorization scales are fixed by the corresponding vector-boson momentum transfer /z^ = 
= -g2 for x = Xi (i = 1, 2). 

Including the SUSY-QCD correction at both qjQjV vertices, the LO order structure 



functions Fi{xj,M'^ 



Fi{xj,M- 



1, . . . , 3 and j = 1, 2) have to be shifted to: 
2a,(/i), 



F,{x,,M' 



1 + 



TT 



^eC{q,,mq,mg) 



(4) 



To illustrate the size of these corrections, we perform a numerical analysis for the light 
scalar Higgs boson h in the decoupling limit of large pseudoscalar masses. Ma ~ 1 TeV. In 
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this case the hght h boson couphngs to standard particles approach the SM values. The only 
relevant processes are then the Higgs-strahlung process qq hV, the vector boson fusion 
mechanism qq qqV*V* —>■ qqh and the gluon-gluon fusion mechanism gg h. 
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Figure 2: Relative corrections due to virtual squark and gluino exchange diagrams to 
Higgs boson production via Higgs-strahlung qq ^ h + W/Z and vector boson fusion qq — > 
qqV*V* qqh [V = W, Z] at the LHC (left) and the Tevatron (right). 



We evaluated the Higgs mass for tan/5 = 30, Ma = 1 TeV and vanishing mixing in 
the stop sector; this yields a value = 112.6 GeV for the light scalar Higgs mass. For 
the sake of simplicity we decompose the K factors K = ctnlo/o'lo into the usual QCD 
part Kqcd and the additional SUSY correction 5susy'- K = Kqcd + ^susy- The NLO (LO) 
cross sections are convoluted with CTEQ4M (CTEQ4L) parton densities and NLO (LO) 
strong couplings a^. The additional SUSY-QCD corrections Ssusy are presented in Fig. || 
as a function of a common squark mass for a fixed gluino mass rrig = 200 GeV [for the sake 
of simplicity we kept the stop mass fixed for the determination of the Higgs mass Mh and 
varied the loop-squark mass independently]. 

The SUSY-QCD corrections increase the Higgs-strahlung cross sections by less than 1.5%, 
while they decrease the vector boson fusion cross section by less than 0.5%. The maximal 
shifts are obtained for small values of the squark masses of about 100 GeV, which are already 
ruled out by present Tevatron analyses [^]; for more reasonable values of these masses, the 
corrections are even smaller. Thus, the additional SUSY-QCD corrections, which are of 
similar size at the LHC and the Tevatron, turn out to be small. For large squark/gluino 
masses they become even smaller due to the decoupling of these particles, as can be inferred 
from the upper squark mass range in Fig. |^. 

In summary, the SUSY-QCD corrections to Higgs boson production in these channels 
are very small and can be safely neglected. 
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3 Associated Higgs production with bb pairs 
3.1 Constraints on the MSSM parameter space 

In the MSSM, the Yukawa couphngs between the Higgs bosons and the down-type fermions, 
in particular the relatively heavy bottom quarks, are enhanced for large tan (3 values. This en- 
hancement can be so significant that it renders the cross section of the associated production 
channel {pp, pp — >■ (p^bb, with (p^ = h^, ^4°, the highest at the Tevatron and the LHC, along 
with the cross section of the gluon fusion mechanism pp gg{Yia, heavy (s)fermion loop) — > 
(fPX 0. The Higgs bosons in this regime decay mainly into hb pairs, leading to 4 6-jets 
which can be tagged experimentally ||29|. Due to the lack of phase space and the reduced 
couplings, the associated production with top quarks is not feasible at the Tevatron, and is 
difficult at the LHC. This makes it possible for the Tevatron Runll and LHC to discover 
Higgs bosons in the (fPbb process and to impose stringent constraints on the SUSY-Higgs 
sector in a relatively model independent way. [At the LHC, the associated H/ A + bb pro- 
duction with the T^T~ and fi~^fj^~ Higgs decay channels is very important |Q, |^ and allows 
to cover most of the parameter space for large tan /?.] 

In Ref. [^, an effective search strategy was presented for the extraction of the signal 



from the backgrounds [which have been calculated]. Using HDECAY [^ to calculate the 
Higgs [and SUSY] spectrum and branching fractions, and combining signals from the search 
of more than one scalar boson [provided their masses differ by less than a resolution Arriexp 
which can be chosen as the total Higgs decay width], contours in the tan/^-m^ plane of the 
MSSM, for which the Tevatron and LHC are sensitive, can be derived. When scanning over 
the parameter space, the set of soft breaking input parameters should be compatible with 
the current data from LEP II and the Tevatron while, preferably, not exceeding 1 TeV. The 
most important parameters here are the masses and mixing of top squarks, and the value 
and sign of the Higgsino mass parameter /z. 

For soft breaking parameters Msoft = = 500 GeV, Fig. ^ shows the 95% C.L. exclusion 
contours in the tan/^-m^ plane, derived from the measurement of a{pp,pp cfPbb bbbb). 
The areas above the four boundaries are accessible at the Tevatron Runll with the indicated 
luminosities and for the LHC with 100 fb"'^. The potential of hadron colliders with these 
processes is compared in Fig. with that of LEPII [where Higgs bosons are searched for in 
the Zh and hA production channels] for the "benchmark" parameter scan "LEPII Scan A2" 
discussed in [^ for ^/s = 200 GeV and a luminosity of 100 pb~^ per experiment. As can 
be seen, the Tevatron can already cover a substantial region with only a 2 fb~^ luminosity. 
Furthermore, for niA ~ 100 GeV, Tevatron and LEPII are complementary. The LHC can 
further probe the MSSM down to values tan/3 ~ 7 (15) for tua < 400 (1000) GeV. 

In conclusion, detecting the (p^bb signal at hadron colliders could effectively probe the 
MSSM Higgs sector, especially for large tan/3 values^. Similar conclusions are reached in 



Ref. ||3^ for the LHC and in Ref. [0, The results given here show a substantial improve- 



ment compared to Ref. [^, where only the pp — >■ (p^bb t^t bb process is discussed at the 

^Note that so far, existing experimental studies are not confirniing the potential of this channel at the 
LHC Q , while the results seem to be more promising at the Tevatron Run II ||] . 
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Figure 3: 95% CL. theoretical estimates of sensitivity contours in the tan/J-m^ planes of 
the MSSM. The areas above the four boundaries can be excluded by the Tevatron Runll and 



the LHC; M,oft = 500 GeV (a) and the "LEP II Scan A2" (b) are shown. From Ref. p| 



Tevatron Runl. Detailed interpretation of the above results in the MSSM and other models 
[such as composite Higgs models with strong dynamics associated with heavy quarks] can be 
found in Ref. |^0[. The analyses can be improved in many ways, for instance with a better 



6-trigger, which bears central significance for the detection of the 6-jets. 



3.2 QCD corrections to Higgs production in heavy quark fusion 

Recently it was proposed that, due to the top-mass enhanced flavor mixing Yukawa coupling 
of the charm and bottom to charged scalar or pseudoscalar bosons (0^), the s-channel 
partonic process c6, cb 0^ can be an important mechanism for the production of 0^ . 
This mechanism is also important for s-channel neutral scalar production via \k fusion^ In 
this section, we describe the complete NLO QCD corrections to these processes. The results 
were originally calculated in Ref. [37|, to which we refer for details. The QCD corrections 



for the SM Higgs production bb ^ H has been also discussed in Ref. [^]. The overlapping 
parts of the two calculations are in agreement. 

The NLO contributions to the process bb — 0° contain three parts: (i) the one-loop 
Yukawa vertex and quark self-energy corrections (Fig. ^o-d); (ii) real gluon emission in qq' 
annihilation (Fig. (iii) s- and t-channel gluon-quark fusion (Fig. ^f-g). In addition, the 
renormalization of the fermion-higgs-fermion Yukawa coupling has to be performed. Since 
the factorization scale /ip = tt^4> is much larger than the mass of the bottom quark, when 
computing the Wilson coefficient functions the 6-quarks were treated as massless partons in 
the proton or anti-proton, similarly to Ref. [^. The only effect of the heavy quark mass 
is to determine at which scale this heavy parton becomes active. (This is the CoUins- 
Wilczek-Zee (CWZ) H scheme). The CTEQ4 PDFs are used to calculate the rates. 



•^Note that the subprocess hh (ffi alone overestnTiates the complete cross section via bottom fusion; one 
has to add consistently the cross sections for hg b4P and gg hhcjp to have a reliable value. 
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because they are consistent with the scheme used in the current study |^ . 




Figure 4: Representative diagrams for charged or neutral (pseudo-) scalar (dashed line) 
production from quark- anti quark and quark-gluon collisions at 0(a;°) and 0{al): (a) leading 
order contribution; (b-d) self-energy and vertex corrections (with counter term); (e) real 
gluon radiation in gg'-fusion; (f-g) s- and t-channel gluon-quark fusion. 

The as corrections involve the contributions from the emission of real gluons, and as 
a result the scalar particle will acquire a non- vanishing transverse momentum Qt- When 
the emitted gluons are soft, they generate large logarithmic contributions of the (lowest 
order) form agin™ {Q^ /Qt)/Qt' where Q is the invariant mass of the scalar and m = 0, 1. 
These large logarithms spoil the convergence of the perturbative series, and falsify the 0{as) 
prediction of the transverse momentum when Qt <C Q- To predict the Qt distribution 
one can use the CoUins-Soper-Sterman (CSS) formalism |^2|, resumming the logarithms of 



the type a" In"^ (Q'^/Qt) /Q'ti orders n in (m = 0, ...,2n — 1). The resummation 

calculation is performed along the same lines as for vector boson production (cf. |E3|). To 
recover the 0{as) cross section, the Wilson coefficients C^^ are included in the resummed 
calculation in The non-perturbative sector of the CSS resummation is assumed to be 
the same as for vector boson production in Ref. |43[| . 

The resummed total rate is the same as the 0{as) rate, when we include C^^ and the 
usual fixed order NLO corrections at high Qt, and switch from the resummed distribution 
to the fixed order one at Qt = Q- When calculating the total rate, we have applied this 
matching prescription. In the case of the scalar production, the matching takes place at 
high Qt ~ Q values, and the above matching prescription is numerically irrelevant when 
calculating the total rate, since the cross sections around Qt ~ Q are negligible. Thus, as 
expected, the resummed total rate differs from the 0{as) rate only by a few percent. Since 
the difference of the resummed and fixed order rates and the /^-factors (cf. Fig. ^) is small, 
we can conclude that for inclusive scalar production once the resummation is performed, the 
0{a'^) corrections are likely to be smaller than the uncertainty from the PDF's. 

Since the QCD corrections are universal, the application to the production of neutral 
scalar or pseudo-scalar 0° via the bb fusion is straightforward. In the following, we will 
consider only the production of the pseudo-scalar within the context of the MSSM. The 
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Figure 5: Cross sections for production in the MSSM with tan/? = 40 at the Tevatron 
and the LHC. (a) The NLO cross sections with the resummed running (sohd) and one-loop 
Yukawa couphng (dashed), as well as the LO cross sections with resummed running (dotted) 
and tree-level Yukawa coupling (dash-dotted) are shown. The cross sections at ^/s = 1.8 TeV 
(thin set of lowest curves) are multiplied by 0.1 not to overlap with the ^/s = 2 TeV curves, 
(b) The NLO (solid), the bb (dashed) and bg (dash-dotted) sub-contributions, and the LO 
(dotted) contributions. The (negative) bg cross sections are multiplied by —1. (c) The NLO 
cross sections with QCD running Yukawa coupling (solid curves) and those with additional 
SUSY corrections (top/bottom dashed lines for yU = +/ — 500 GeV). 

total LO and NLO cross sections for the inclusive processes pp,pp —>■ A^X at the Tevatron 
and the LHC are shown in Fig. |^ for tan f3 = 40. For other values the cross sections can be 
obtained by scaling with the factor (tan/?/40)^. 

Fig. ^ shows a significant improvement from the pure LO results (dash-dotted curves) 
due to the resummation of the large logarithms of m^/ml into the running coupling. The 
good agreement between the LO results with running coupling and the NLO results is due to a 
non-trivial, and pro cess- dependent, cancellation between the individual 0{as) contributions 
of the bb and bg sub-processes (which are connected via mass factorization). 

For large tan/3, the SUSY correction to the running cfp-b-b Yukawa coupling can be 
significant |^4[, and can be included in a similar way as it is done for the (fPbb associate 
production [|^. To illustrate the effects of these corrections, all MSSM soft-breaking pa- 
rameters and /i were set to 500 GeV. Depending on the sign of /x, the correction to the 



coupling can take either the same or opposite sign as the full NLO QCD correction [45 
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Figure 6: The /^-factors for production in the MSSM with tan/5 = 40 for the NLO 
{K = (Tnlo/o-lo, sohd hues), bb {K = cr;,5/aLo = (c^lo + '^<7bb)/<^LO, dashed hues), and bg 
{K = —Aabg/a-Lo, dash-dotted hues) contributions, at the Tevatron (a) and LHC (b). 

In Fig. I^c, the sohd curves represent the NLO cross sections with QCD correction alone, 
while the results including the SUSY corrections to the running bottom Yukawa coupling are 
shown for /x = +500 GeV (top dashed curves) and fi = —500 GeV (bottom dashed curves). 
These partial SUSY corrections can change the cross sections by about a factor of 2. 

The i^'-factors, the ratios of the NLO versus LO cross sections as defined in Ref. [ p7[] , 
for the pp,pp A^X processes are presented in Fig. ^ for the MSSM with tan/3 = 40. 
Depending on the A^ mass, they range from about — (16~17)% to +5% at the Tevatron and 
the LHC. The uncertainties of the CTEQ4 PDFs for A°-production at the Tevatron and the 
LHC are summarized in Fig. 0. 

The transverse momentum distributions of A^, produced at the upgraded Tevatron and 
at the LHC, are shown in Fig. ^ for various A'^ masses with tan /3 = 40. The solid curves are 
the result of the multiple soft-gluon resummation, and the dashed ones are from the 0{as) 
calculation. The fixed order distributions are singular as Qt — >■ 0, while the resummed ones 
have a maximum at some finite Qt, and vanish at Qt = 0. When Qt becomes large, of the 
order of rriA, the resummed curves merge into the fixed order ones. The average resummed 
Qt varies between 25 and 30 (40 and 60) GeV in the 200 to 300 (250 to 550) GeV mass 
range of rriA at the Tevatron (LHC). 

In summary, the overall NLO corrections to the pp,pp A^X processes are found to 
vary between — (16~17)% and +5% at the Tevatron and the LHC in the relevant range of 
the A^ mass. The uncertainties of the NLO rates due to the different PDFs also have been 
systematically examined, and found to be around 20%. The QCD resummation, including 
the effects of multiple soft-gluon radiation, was also performed to provide a better prediction 
of the transverse momentum distribution of the scalar 0°. This latter is important when 
extracting the experimental signals. Similar results can be easily obtained for the other 
neutral higgs bosons {h^ and H^) by properly rescaling the coupling. These QCD corrections 



29 



= 1.2 

O 1 ,1 5 

to 

1 .1 



Tevatron 
VS = 2 TeV 



CTEQ4HJ ..' 



CTEQ4-M 



CTEQ4A5 



CTEQ4A1 



(o) 



LHC 
VS = 1 4 TeV 



CTEQ4A5 



t2< 



CTEQ4A1 



(b) - 



200 400 600 SOO 1 OOO 

(GeV) 



Figure 7: The ratios of NLO cross sections computed by four different sets of CTEQ4 PDFs 
to the cross section computed by CTEQ4M for neutral pseudo-scalar {A^) production in the 
MSSM with tan/5 = 40, at the upgraded Tevatron (a) and the LHC (b). 
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Figure 8: Transverse momentum distributions of pseudo-scalar produced via hadron 
collisions, calculated in the MSSM with tan/5 = 40. The resummed (solid) and 0{as) 
(dashed) curves are shown for = 200, 250, and 300 GeV at the upgraded Tevatron (a), 
and for tua = 250, 400, and 550 GeV at the LHC (b). 



can also be applied to the generic two higgs-doublet model (called type-III 2HDM[ 
which the two higgs doublets $i and $2 couple to both up- and down-type quarks. 
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4 Higgs search in the 77+jet channel at LHC 



The observation of a Higgs boson with a mass Mh < 140 GeV at the LHC in the inclusive 
channel — 77 + X is not easy |57, Q as it is necessary to separate a rather elusive Higgs 
boson signal from the continuum background. In Ref. |^ the reaction pp H{-^ 77)+jet, 
when the Higgs boson is produced with large transverse momentum recoiling against a hard 
jet, was analyzed as a discovery channel. The signal rate is much smaller, but there remains 
enough events to discover the Higgs boson at a low luminosity LHC. It is important to 
note that the situation with the background is undoubtedly much better in the case of Higgs 
production at highpT-. Thus, one has S/B 1/2 — 1/3 for CMS and ATLAS correspondingly, 
providing a discovery significance of 5 already with an integrated luminosity of 30 fb~^. 
Furthermore, recent achievements in calculations of QCD next-to-leading corrections have 
shown an enhancement of the signal against the background. This circumstance together 
with the possibility to exploit the event kinematics in a more efficient way allow the hope 
that this reaction will be the most reliable discovery channel for Higgs bosons with masses 
Mh = 110- 135 GeV. 

Typical acceptances of the LHC detectors ATLAS and CMS were taken into account in 
the analysis: two photons are required with p] > 40 GeV for each photon (harder than for 
the inclusive channel), and \ri^\ < 2.5, while a jet was required with Ej'^* > 30 GeV and 
|?7jei| < 4.5, thus involving the forward parts of the hadronic calorimeter. The isolation cut 
AR > 0.3 was applied for each 77 and 'yg{q) pair. 

There are three QCD subprocesses giving a signal from the Higgs boson in the channel 
under discussion in QCD leading order: gg ^ H + g, gq H + q and qq ^ H + g. It was 
found that the gg^H + g^'j'y + g subprocess gives the main contribution to the signal 
rate. In total, the QCD signal subprocesses give 5.5, 10.6 and 9.8 fb for Mh = 100, 120 and 
140 GeV, correspondingly within the kinematical cuts described above. 

Another group of signal subprocesses includes the electroweak reactions of Higgs produc- 
tion through WW or ZZ fusion and in association with W or Z boson, where one should 
veto the second quark jet. The EW signal rate is at the level of 10% of the QCD signal. 

Both the reducible and irreducible backgrounds, pp 77+jet have been discussed in the 
QCD section of these Proceedings. It was found that in total it is about 19, 31 and 32 fb in 
the 1 GeV bin for Mh = 100, 120 and 140 GeV, correspondingly. 

Further improvement of the S/B ratio can be obtained by studying the kinematical 
distributions of the 3-body final states in the subprocesses under discussion. The background 
processes contribute at a smaller a/^ in comparison to the QCD signal processes. So, the 
corresponding cut improves the S/B ratio: e.g., the cut > 300 GeV suppresses the 
background by a factor of 8.7 while the QCD signal is suppressed only by a factor of 2.6. 
This effect is connected with the different shapes. Fig. |^, of the jet angular distributions in the 
partonic c.m.s. for the signal and background. Indeed, for the dominant signal subprocess 
gg ^ H + g, a. set of possible in spin states does not include spin 1, while the spin of the 
out state is determined by the gluon. It means, in particular, that the S-wave does not 
contribute here. At the same time, in the dominant background subprocesses gq —>■ 'J'J + q 
and qq ^ 'J'J + g, the same spin configurations are possible for both in and out states. It 
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was found that the cut on the partonic colhsion energy matches this spin-states effect, 
and the best S/ B ratio is obtained at > 300 GeV. One can try to exploit this effect to 
enhance the signal significance with the same level of the S/B ratio. Indeed, if one applies 
the cut on the angle between the jet and the photon in partonic c.m.s. cos'(9*_^ < —0.87 for 
< 300 GeV and add such events to the events respecting the only cut > 300 GeV, 
then the S/ B change is rather small, while the significance is improved by a factor of about 
1.3. The same effect can be observed with the cut on the jet production angle in the partonic 



c.m.s. "i^jg^, but one should note that the two variables. 



e 



and '(9*gj, are correlated. 



It is 



desirable to perform a multivariable optimization of the event selection. 

Note that this is a result of a LO analysis, the task for the next step is to understand 
how this effect will work in presence of NLO corrections to both the signal and background. 



In the analysis performed in Ref. [59| the factor K^^^ = 1.6 was used to take into 
account the QCD next-to-leading corrections for both the signal and background subpro- 
cesses. In Ref. [^0|, |51|, this assumption was confirmed by an accurate evaluation of NLO 
corrections to the signal subprocesses (where for the evaluation of the two-loop diagrams, 
the effective point-like vertices were used in the limit Mh -C rrit |^). For the background, 
the corresponding analysis has shown that the NLO corrections are not larger than 
50%. Thus, an attractive feature of the pp H{-^ 77)+jet channel is that theoretical 
uncertainties related to higher order QCD corrections can be under control. 



fb/deg 




120 140 160 ISO 




20 40 60 80 100 120 140 160 




20 40 60 80 100 120 140 160 180 





20 40 60 80 100 120 140 160 180 



31 



Figure 9: Distributions in the jet production angle 'd*^^^ and the angle 'd*^ between jet and 
the photon with smaller pt in partonic c.m.s. for the QCD signal (S) and background (B). 
The Higgs mass is taken to be = 120 GeV. Upper plots - no cut, in others > 210 
and 300 GeV correspondingly. The M^^ bin is equal to 1 GeV. 
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Signatures of Heavy Charged Higgs Bosons at the LHC 

K.A. ASSAMAGAN, A. DjOUADI, M. Drees, M. Guchait, R. Kinnunen 
J.L. Kneur, D.J. Miller, S. Moretti, K. Odagiri and D.P. Roy 

Abstract 

We analyze the signatures of the charged Higgs particles of the Minimal Supersym- 
metric extension of the Standard Model at the LHC. We will mainly focus on the large 
Mff± range where the charged Higgs boson is produced through the gluon-bottom or 
gluon-gluon mechanisms. The resulting signal is analyzed in its dominant — > tb 
as well as subdominant decay channels. Simulations for the detection of the charged 
Higgs boson signals in the decay channels — > t^u and — > cs,W^h or tb are 
performed in the framework of the CMS and ATLAS detectors, respectively. 

1 Introduction 

The minimal supersymmetric Standard Model (MSSM) contains two complex Higgs doublets, 
01 and 02, corresponding to eight scalar states. Three of these are absorbed as Goldstone 
bosons leaving five physical states - the two neutral scalars H^), a pseudo-scalar (A^) and 
a pair of charged Higgs bosons (H^). All the tree-level masses and couplings of these particles 
are given in terms of two parameters, and tan/5, the latter representing the ratio of 

the two vacuum expectation values While any one of the above neutral Higgs bosons 
may be hard to distinguish from that of the Standard Model, the carries a distinctive 
hall-mark of the SUSY Higgs sector. Moreover the couplings of the are uniquely related 
to tan P, since the physical charged Higgs boson corresponds to the combination 

= -(f)fsm(3 + (j)fcos(3. (1) 

Therefore the detection of H"^ and measurement of its mass and couplings are expected to 
play a very important role in probing the SUSY Higgs sector. 

The search for charged Higgs bosons is one of the major tasks of present and future high- 
energy colliders. In a model independent way, LEP2 has set a lower limit on the mass, 
'^Hi ^74 GeV, for any value of tan/? [Q. At the Tevatron, the CDF and D0 collaborations 
searched for bosons in top decays through the process pp ti, with at least one of the 
top quarks decaying via t H^h, leading to a surplus of r's due to the tv decay; 

they excluded the low and high tan/? regions [where the branching ratios for this decay is 
large] almost up to the Mii± ~ nit limit 0. Detailed analyses at the LHC have shown that 
the entire range of tan /5 values should be covered for M//± < rrit using this process. 

At this workshop, we focused on the large mass region, Mii± > rrit, where the previous 
production process is not at work and for which only a few preliminary studies have been 
performed. We summarize our work in this contribution. After a brief summary of the 
decay modes [both in the MSSM and in some of its extensions], we will discuss in section 
3, the various signals for a heavy charged Higgs boson at the LHC. We will then present, in 
sections 4 and 5, two simulations for the detection of the signals in the decay channels 
T^u in the CMS and cs, W^h, tb in the ATLAS detectors. 
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2 Production and decay modes of the bosons 

The decays of the charged Higgs bosons are in general controlled by their Yukawa couplings 
to up- and down-type fermions u, d given by 

U+ [cot (3 rriu UidjL + tan (3 nid Uidj^] , (2) 



w 



For values tan/3 > 1, as is the case in the MSSM, the couplings to down-type fermions are 
enhanced. The coupling H~tb, which is of utmost importance in the production and the 
decays^ of the bosons, is large for tan (3 ~ 1 and ~ mt/mh. Interestingly these two 
regions of tan/3 are favored by 6-r unification for a related reason: i.e. one needs a large 
tbH"^ Yukawa coupling contribution to the RGE to control the rise of nib as one goes down 
from the GUT to the low energy scale M. 
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Figure 1: Branching ratios of the charged Higgs boson decays for tan/3 
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2 and 30. They 



are obtained using the program HDECAY |10 . 



The value of tan /3 determines to a large extent the decay pattern of the charged Higgs 
bosons. For large tan /3 values the pattern is simple, a result of the strong enhancement of the 
couplings to down-type fermions: below the top-bottom threshold, bosons will decay 
into TUr pairs while above this threshold, they will decay into tb pairs with BR ~ 85% and 

"^It should be mentioned that most analyses of the boson decay modes and detection signals at 
colliders are based on the lowest order vertex, represented by the Yukawa coupling of eq. (|^), but improved 
by standard QCD corrections ||^ by using the running quark masses. One loop electroweak corrections to 
this vertex can give a large variation in the signal cross-section at high or low tan/3, as recently shown 
in p. The corresponding correction from SUSY-QCD loops is possibly large j^] depending on the SUSY 
parameters [but for the production, they are not yet completely available] . The inclusion of these corrections 
is evidently important for a quantitative evaluation of the signal. 
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TVr pairs with BR ~ 15% for large enough Mii± values. For small tan/5 values, tan/? < 5, 
the pattern is more complicated, in particular around and below the th threshold. Decays 
into Wh final states play an important role since they reach the level of several ten percent 
leading to a significant reduction of the dominant branching ratio into tv states. Note that 
the off-shell three body decays ^ ^ bt* -> bbW^ and hW\AW* hff,Aff 

[the latter being kinematically forbidden at the two-body level] can be rather important. 
The branching ratios are summarized in Fig. 1 for the values tan/5 = 2 and 30. 



In the MSSM, the charged and pseudoscalar Higgs boson masses are related 



M^± = Ml + M, 



w 



(3) 



and the LEP limit on the lightest scalar and pseudoscalar Higgs masses, mho^niAo) ^ 90 — 
100 GeV implies first, that Mh± > 120 GeV [Mh± > 200 GeV for tan/3 = 2] and second, 
that the — »• Wh^{WA^) decay channel has as high a threshold as the tb channel, while 
the latter has a more favorable coupling. Consequently the —>■ Wh^{WA^) decay BR 
is restricted to be < 5% over the LEP allowed region [Fig. 1]. However the constraints 
discussed above do not hold in singlet extensions of the MSSM like the NMSSM |11]|. Con- 
sequently — > Wh^{WA^) can be the dominant decay mode for Mh± ~ 160 GeV in the 
low tan/3 region and lead to a spectacular signal at the LHC, as illustrated in Table 1. This 
decay channel will be analyzed in detail in the next sections. 

Table 1. Maximal branching fractions for W{h^, A^) decay in the NMSSM for fixed input 

values of tan f3 and output mass of ~ 160 GeV. The values of the /i?, ^? masses and branching 
fractions are shown along with the corresponding model parameters. Also shown are the t — > bH^ 
branching fraction and the size of the resulting W{hi, A^) decay signal at LHC. 
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An important point which should be mentioned, is that in most of the analyses of the 
H"^ signals, it is always assumed that it decays only into standard particles and that the 
SUSY decay modes are shut. But for large values of Mh±, at least the decays into the 
lightest neutralinos and charginos [and possibly into light sleptons and t, b squarks] can be 
kinematically allowed. These modes could have large decays widths, and could thus suppress 
the — »• tb branching ratio in a drastic way 

In Fig. 2, the branching fraction BR(if^ — > x'iXf ) [with i =1-4 and j=l-2] are shown 
as function of Mfj± for the four values tan/5 = 2,5, 10 and 30. The choice of the gaugino 
and higgsino mass parameters M2 = /i = 200 GeV has been made leading to the lightest 
chargino and neutralino masses m^o ~ 80-90 GeV and ~ 125-150 GeV depending on 
the value of tan f3 [small masses are obtained with small tan f] input] . The values of the scalar 
masses are such that sleptons and squarks are too heavy to appear in the decay products of 
the boson. As can be seen, for small and large values of tan (3, the H^tb couplings are 
enhanced and the chargino/neutralino decays are important only for large masses where 
many channels are open. For intermediate tan/5 values, the H^tb Yukawa couplings 
are suppressed, and the chargino/neutralino decays are dominant for charged Higgs boson 
masses of a few hundred GeV. 

In scenarii where sleptons and squarks [in particular stop and sbottom squarks] are 
also light, if^ bosons decays into these states might be kinematically possible as well and 
would be dominant. This will again suppress in a dramatic way the branching ratio for the 
if^ — > tb signature |T^. These SUSY decays, although discussed in the literature, have not 
been analyzed experimentally up to now. They should, however, not be overlooked for heavy 
charged Higgs bosons, as they might jeopardize the detection of these particles at the LHC. 
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Figure 2: Branching ratios of the charged Higgs boson decays into charginos and neutralinos 
as a function of Mfj± for a set of tan /3 values; M2 and /i are fixed to 200 GeV. 
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Finally, we briefly discuss the production modes of a heavy charged Higgs boson, with 
mfj± > rrit, at the LHC The two mechanisms which have sizeable cross sections are: 



pp 
pp 



99 1 m' - 



tH- {tH+) 
tH-b + tH+b 



13, |14 
B 0, 



(4) 



The signal cross-section from the 2^2 mechanism gb — >■ tH~ [where the b quark is obtained 
from the proton] is 2-3 times larger than the 2 — > 3 process 99 /qq — tbH~ [where the H~ 
boson is radiated from a heavy quark line]. This is shown in Fig. 3a at LHC energies for 
the values tan/3 = 2 and 40. When the decays H'^ —>■ tb and t — > Wb take place, the first 
process gives rise to 3 6-quarks in the final state while the second one gives 4 6-quarks. Both 
processes contribute to the inclusive production where at most 3 final 6-quarks are tagged. 



However, the two processes have to be properly combined [|T^ to avoid double counting of 
the contribution where a gluon gives rise to a hb pair that is collinear to the initial proton. 
The cross section of the inclusive process in this case is shown in Fig. 3b, and is mid-way 
between the two cross sections eqs. (4) 
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Figure 3: Production cross sections for charged Higgs bosons at the LHC for tan [3 = 2 and 
40. (a) Individual cross sections from the 99 /qq and gb processes and (b) combination of the 
two processes with the subtraction of the common piece. 

Other mechanisms for production at hadron colliders are the Drell-Yan type process 
for pair production, qq H^H^ , the associated production process with W bosons, qq 
H^W"^ |TP[ and the gluon-gluon fusion process for pair production, gg — > H~^H~~ [^. 
However, the rates are rather small at the LHC, in the first case because of the weak couplings 
and the low quark luminosities at high energies and in the second case because the process is 
induced by loops of heavy quarks and is thus suppressed by additional electroweak coupling 
factors. We will thus focus in this study on the two processes eq. (4). 
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3 Signatures of bosons at the LHC 

The t —>■ hH^ decay is known to provide a promising signature for charged Higgs boson search 
at the LHC for < rrit. But it is hard to extend the search beyond rrit, because in this 
case the combination of dominant production and decay channel, tH~ —>■ ttb, suffers from a 
large QCD background 15]. Moreover the subdominant production channels of H^W^ 
and H^H^ have been found to give no viable signature at LHC [0. In view of this we 
have undertaken a systematic study of a heavy signature at the LHC from its dominant 
production channel gb{gb) tH~(tH'^), followed by the decays H~ — >• ib,Tu and W~h^. 
While the first decay represents the dominant channel of charged Higgs bosons, the tu and 
Wh'^ are the largest subdominant channels in the high and low tan /? regions respectively, 
with [see also Fig. 1] 

(tan/? > 10) ~ 15% and Bwhoitan/3 = 1 - 5) < 5%. (5) 

The signature for the dominant decay channel of H~ —>■ tb has been analyzed separately 
assuming three and four 6-jet tags. The analyses presented in this section are based on 
parton level Monte Carlo programs with a Gaussian smearing of lepton and jet momenta for 
simulating the detector resolution. 



(i) H — > tb Signature with Four 6-tags []T 



The dominant signal and background processes are 

gg tH'b + h.c. ttbb, (6) 

gg ^ ttbb, (7) 
followed by the leptonic decay of one top and hadronic decay of the other, i.e. 

ttbb bbbbW+W- bbbUvqq. (8) 

A basic set of kinematic and isolation cuts, 

Pt > 20 GeV, \ri\ < 2.5, AR = [{A^f + {Arif]^^^ > 0.4 (9) 

is imposed on all the jet and lepton momenta. The pt cut is also imposed on the missing-p^-, 
obtained by vector addition of the p^'s after resolution smearing. This is followed by the 
mass reconstruction of the W and the top quark pair, so that one can identify the pair 
of 6-jets accompanying the latter. While the harder of these two 6-jets comes from 
decay in the signal, both of them come mainly from gluon splitting in the background. 
Consequently the S/B ratio is improved by imposing the following cuts on this 6-jet pair: 

Mbb > 120 GeV, Eb, > 120 GeV and cos^;,;, < 0.75. (10) 



^While this work was initiated earlier, some of the issues analyzed during the workshop led to the final 
version presented here. 
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Then each of this 6-jet pair is combined with each of the reconstructed pair of top to give 4 
entries for the invariant mass Mj^ per event. One of these 4 entries corresponds to the if^ 
mass for the signal event, while the others constitute a combinatorial background. Fig. 4 
shows this tb invariant mass distribution for the signal @ and background (|^). The right 
hand scale corresponds to the cross-section for = 0.1 - i.e. an optimistic 6-tagging efficiency 
of efc = 0.56. Reducing it to a more conservative value of = 0.4 would reduce both the 
signal and background by a factor of 4 each. 

b t H + c.c. ^bbtt^4b2jl + missing energy/momentum 

0.30 
0.25 



^ 0.20 

\ 

* 

t 0.15 



§ 0.10 

\ 
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0.05 



0.00 

200 400 600 800 1000 

iliin) (GeV) 

Fig. 4: The reconstructed tb invariant mass distribution of the signal (^) and the QCD 
background (|^) in the isolated lepton plus multi-jet channel with 4 6-tags. The scale 
on the right corresponds to a 6-tagging efficiency factor ef = 0.1. 

Table 2. Number of signal and background events in the 4 6-tagged channel per 100 fb^^ luminosity 
in a mass window of Mjj± it 40 GeV at tan /? = 40 (e;, = 0.4). 
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Table 2 lists the number of signal and background events for a typical annual luminosity 
of 100 fb~^, expected from the high luminosity LHC run, assuming eb = 0.4. While the S/B 





t 

\ 


MRS-L0(05A) 


[ Q = M 
Mh+ = 


= -JS] 

310 GeV 






ii 
il 

1 




M„+ = 
Mh+ = 
Mh+ = 


407 GeV 
506 GeV 
605 GeV 




' 1 j- 


1 
1 

1 

^ 1 

^ 1 




Mh+ = 
Mh+ = 


704 GeV 
803 GeV 






^ 1 




tan|3 = 


40. 




' 1 / 


\ 


1 





























0.020 



42 



ratio is > 1, the viabihty of the signal is hmited by the signal size[|. One expects a > Scr 
signal up to Mij± = 600 GeV at tan/3 = 40. The signal size is very similar at tan/5 = 1.5, 
but smaller in between [the signal process is controlled by the tbH^ Yukawa coupling, 
eq. (H), which is large for tan/5 ~ 1 and ~ nit/ nib, as discussed previously]. 



(ii) H"^ — i> tb Signature with Three 6-tags [j2T 



The contributions to this signal come from (|^) as well as 

gb ^ tH' + h.c. ^ tib + h.c, (11) 

followed by the leptonic decay of one top and hadronic decay of the other. The signal cross- 
section from ([TTD is 2-3 times larger than from (|^) [Fig. 3] , while their kinematic distributions 
are very similar. Combining the two cross sections and subtracting the overlapping piece to 
avoid double counting results in a signal cross-section, which is mid-way between the two; 
see Fig. 3. 

The background comes from (|^) as well as 

gb ttb + h.c. and gg tig, (12) 

where the gluon jet in the last case is mis-tagged as a 6-jet. Assuming the standard mis- 
tagging factor of 1% this contribution turns out in fact to be the largest source of the 
background, as we see below. 

The basic kinematic cuts are as in (|^) except for a harder p^-cut, 

Pt > 30 GeV, (13) 

since the 3 6-jets coming from and tt decays are all reasonably hard. This is followed 
by the mass reconstruction of the top quark pair as before, so that one can identify the 
accompanying (3rd) 6-jet. We impose a 

Pt > 80 GeV (14) 

cut on this 6-jet to improve the S/B ratio. Finally this 6-jet is combined with each of 
the reconstructed top pair to give two entries of Mth per event. One of them corresponds 
to the mass for the signal while the other constitutes the combinatorial background. 
Fig. 5 shows this tb invariant mass distribution of the signal along with the above mentioned 
backgrounds, including a 6-tagging efficiency factor of 

eb = 0.4. (15) 

While the S/B ratio is < 1 the signal cross-section is much larger than the previous case. 
Table 3 lists the number of signal and background events for a luminosity of 100 fb^^ at 
tan/? = 40. The results are very similar at tan/5 = 1.5. Comparing this with Table 2 we see 
that the 5*/ \^ ratio is very similar in the two channels. One should bear in mind however 
the larger pt cut (|13D assumed for the 3 6-tagged channel. The cross- sect ions in both the 



cases were calculated with the MRS-LO structure functions 122 



^Increasing the px cut of &-jets from 20 to 30 GeV would reduce the signal (background) size by a factor 
of about 3(4), hence reducing the viability of this signal. 
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Fig. 5: The reconstructed tb invariant mass distribution of the H"^ signal and different QCD 
backgrounds in the isolated lepton plus multi-jet channel with 3 6-tags. 



Table 3. Number of signal and background events in the 3 6-tagged channel per 100 fb luminosity 
in a mass window of Mfj± ± 40 GeV at tan/3 = 40 (e?, = 0.4). 



Mh± (GeV) 


S 


B 


Sy/B 


310 


133 


443 


6.2 


407 


111 


403 


5.6 


506 


73 


266 


4.5 


605 


43 


156 


3.4 



(iii) Tu Signature USsj : 

Following the analysis of Ref. |^ a more exact simulation of a heavy signature in 
the Tu decay channel was done for the CMS detector using PYTHIA [Q. The results will 
be presented in the next section. By exploiting the distinctive r polarization one can get at 
least as good a signature here as in the tb channel for the large tan f3 region. 



(iv) W^h^ Signature [|25 



For simplicity we have estimated the signal cross-section from 

gb ^ tH- + h.c. ^bW+W-h^ + h.c., (16) 
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followed by — » bb, — > iu and — > qq. Thus the final state consists of the same 
particles as the dominant decay mode of eq. (0). Thus we have to consider the background 
from the H~ —>■ tb decay ([TT|) along with those from the QCD processes of eq. ([T^). 

We require 3 6-tags along with the same basic cuts as in section (ii). This is followed by 
the mass reconstruction of and the top, which helps to identify the accompanying 6-pair 
and the W. The resulting bb and Wb invariant masses are then subjected to the constraints, 

Mbb = rriho ± 10 GeV and mwb rrit ±20 GeV. (17) 

The mass constraint and the veto on the second top helps to separate the —>■ W^h^ 
signal from the backgrounds. However the former is severely constrained by the signal size 
as well as the S/ B ratio. Consequently one expects at best a marginal signal in this channel 
and only in a narrow strip of the Mj:f±-tan j3 parameter space, at the boundary of the LEP 
exclusion region. Fig. 6 shows the signal (|T6|) along with the backgrounds from (|1T]) and 
(12) against the reconstructed mass at one such point - Mh± = 220 GeV and tan (3 = 2. 
Note that, as discussed in section 2, in extensions of the MSSM, the H"^ — Wh^(WA^) can 
be the dominant decay mode for Mff± ~ 160 GeV in the low tan/3 region and lead to a 
spectacular signal at the LHC; see Table 1. 



0.0025 



0.002 - 




150 200 250 300 350 400 450 

I^Wbb 

Fig. 6 The H"^ Wh^ signal cross-section is shown against the reconstructed H mass for 
Mh± = 220 GeV and tan (3 = 2 along with the tb and the QCD backgrounds. 

It should be mentioned here that these parton level Monte Carlo analyses of the 
signature in tb and Wh^ decay channels need to be followed up by detailed simulation 
with PYTHIA, including detector acceptance, as in the case of the ru channel discussed 
in the next section. Some work has started here along this line for the ATLAS detector; 
this is summarized in section 5. One should also bear in mind the possibility of large 
radiative corrections to the Yukawa coupling eq. (@); it is evidently important to include 
these corrections for a quantitative evaluation of this signal. 
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4 The H+ TP mode in CMS 



4.1 Introduction 



As mentioned in the previous section, the hadronic r signature of a heavy charged Higgs 
boson from pp — ^ tiJ^ at the LHC is usefuL In this contribution, we study the search of 
heavy H"^ bosons in the CMS detector with a reahstic simulation using the procedure of 
Ref. p3|| to select the events and to exploit the r polarization effects. The main backgrounds 



are due to tt and py+jet events. The VT+jet background can be effectively reduced with W 
and top mass reconstruction and 6-tagging. Although for tt and l^+jet events the transverse 
mass reconstructed from the T-jet and the missing transverse energy is bounded from above 
by the W mass, some leaking of the backgrounds into the signal region can be expected due 
to the experimental resolution of the E^^^^ measurement. 

4.2 Event selection and expectations for CMS 

Events are generated with PYTHIA using the process bg tH^. The results from 
Ref. 1^ with a subtraction of double counting between the gb — > tH^ and gg — >■ tbH^ 
processes are used to normalize the PYTHIA cross sections. The ru branching ratio 



is calculated with the HDECAY program []TU| and used in the simulation. A heavy SUSY 



particle spectrum (1 TeV) is assumed with no stop mixing. The decay matrix elements 



with polarization effects I^S], ^ are added in PYTHIA. For mH± = 400 GeV and tan/3 
= 40 about 1700 signal events, including only one-prong hadronic r decays, are expected 
for an integrated luminosity of 30 fb~^. The jets and the missing transverse energy are 



reconstructed with a fast simulation package CMS JET . For b-tagging, results obtained 



from a full simulation and reconstruction of the CMS tracker are used ||28[| . 

The real r jet is chosen as the r jet candidate requiring > 100 GeV and \rj\ <2.5. 
The events can be triggered with a multi-jet trigger and a higher level r trigger even in 
the high luminosity running conditions. The r selection is performed here using only the 



tracker information. The algorithm of ref. [23| to remove the transverse components of the r 



polarization is used requiring i = p'^ / E'^^^^ > 0.8, where p'^ is the momentum of a hard pion 
from T decay in a cone of AR < 0.1 around the calorimeter jet axis and E'^^'^^ is the hadronic 
energy of the r jet {E^ > 100 GeV) reconstructed in the calorimeters (electromagnetic and 
hadronic) in a cone of AR < 0.4. The efficiency of this r selection for the signal events is 
20% while for the tt events the efficiency is only 0.4% (including the Et threshold for jet). 
A reconstruction efficiency of 95% is assumed for the hard isolated track from r. 

A large missing transverse energy is expected in the signal events due to the neutrino from 
if^ decay. The is reconstructed with the CMSJET package, where the calorimeter 

response is parametrized including the effects of the detector cracks and the volumes of 
degraded response. Efficiency of the cut E™^'^ > 100 GeV is about 75% for the signal events 
and about 39% for the tt background. 

A visible signal for the Higgs can be obtained in the transverse mass reconstructed from 
the T-jet and the missing transverse energy if the hadronic decay of the associated top quark 
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is selected. For the reconstruction of the W and top masses the events with at least three 
jets with Ef > 20 GeV, in addition to the r jet, are selected. The W and top masses are 
reconstructed minimizing the variable = i^jj ~ + (^iii ~ ^tY-, where mw and 

mt are the nominal W and top masses. A Gaussian resolution of 13.6 GeV is found for 
the reconstructed top mass. The fraction of events where the three jets are found and the 
reconstructed W mass is within m^^i 15 GeV and the reconstructed top mass within m^i 
20 GeV is 54% for the signal, 59% for the tt background and 8% for the ly+jet events. 

After the W and top mass reconstruction and the mass window cuts b-tagging is applied 
on the jet not assigned to the W . This jet is required to be harder with El^^ > 30 GeV. The 
tagging efficiencies based on the impact parameter method obtained from a full simulation 
and track reconstruction in the CMS tracker are used p8|. At least two tracks with pt > 
1 GeV and impact parameter significance cr*^ > 2 are required inside the jet reconstruction 
cone of 0.4. For b-jets with Et = 50 GeV the efficiency is found to be ~ 50% averaged over 
the full T] range (|?7| < 2.5). The mis-tagging rate for the corresponding light quark and 
gluon jets is 1.3%. 



70 
60 
50 
40 
30 
20 
10 - 




pp tH^'^.H^ TV, t 



qqb 



m„ = 400 GeV, tan/J = 40 
a) 



Total background 



Signal 



50 100 150 200 250 300 350 400 450 500 
m,(T jet, E,™") (GeV) 



> 








o 




o 


16 


\ 


JD 


14 


Q. 




O 




* 


12 






o 




CO 


10 


c 




<U 




Ev 


8 




6 
4 




2 








pp tH*,H* 



TV, t — > qqb 
m„ = 400 GeV, tanjS = 40 
jet and top veto 
b) 



Total background 



Signa 



300 400 500 

m,(T jet, E,™") (GeV) 



Figure 7: a) Transverse mass reconstructed from r jet and EJ^^'"^^ for H"^ —>■ tu from pp 
tH^ with rriHi = 400 GeV and tan /3 = 40 over the total background from tt and W^+jet 
events, b) the same as in a) but with a veto on a central jet and a second top. 



The reconstructed transverse mass m^^*^ over the total background is shown in Fig. 7a for 
mH± = 400 GeV and tan/3 = 40 for 30 fb"^ For m^" >100 GeV about 44 signal events 
are expected for m^^i = 400 GeV and tan/5 = 40 and about 25 events for mH± = 200 GeV 
and tan/3 = 30, for an integrated luminosity of 30 fb"^. About 5 background events from 
tt and W + jet are expected for mJji' >100 GeV. Further reduction of the ti background is 
still possible using a jet veto cut and a veto on a second top in the event. Since a soft and 
a relatively forward spectator b-jet from the production process is expected in the signal 
events, a central and hard jet veto with \ri^^^\ < 2 and i?^^* > 50 GeV is used. For the 
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reconstruction of the second top from the t jet, missing energy and one of the remaining 
jets, the longitudinal component of the missing energy is first resolved from the W mass 
constraint selecting the smaller of the two solutions. The reconstructed top mass is required 
to fall outside the window of mt± 60 GeV. The central jet veto and the second top veto, 
being closely correlated cuts, reduce ti background by a factor of ~7. The efficiency for the 
signal is 54%. The transverse mass distribution over the total background including 
the jet and second top veto is shown in Fig. 7b for mH± = 400 GeV and tan (3 = 40 and in 
Fig. 8a for m^i = 200 GeV and tan f3 = 30. 
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Figure 8: a) Transverse mass reconstructed from r jet and E^^^^ for — ^ tu from pp — >■ 
tH^ with mH± = 200 GeV and tan /? = 30 over the total background from tt and W^+jet 
events with central jet and second top veto, b) the same as in a) but for A0 > 50° where 
A0 is the angle between the r jet and the vector in the transverse plane. 

The visibility of the signal can be significantly improved, especially at m//± = 200 GeV, 
with a cut on the A0 angle between the r jet and the EY^^^" . Although A0 is directly 
proportional to m^-^, a cut in A0 suppresses the background efficiently at the lower end of 
the expected signal region as can be seen from Fig. 8b showing the signal over the total 
background with A0 > 50° for mH±= 200 GeV and tan/3 = 30. 



4.3 Conclusion 

Our preliminary study leads to the conclusion that tv from pp —>■ tH^ is a promis- 

ing discovery channel for charged Higgs bosons at the LHC. For the evaluation of the final 
discovery reach in the mA, tan (3 parameter space a detailed simulation of the measure- 
ment for the background events is needed. The study can be extended to high luminosity 
but some additional loss of efficiency should be expected due to the harder Ej^^ cuts due to 
trigger requirements. 
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5 The H+ cs, Wh, tb modes in ATLAS 
5.1 Introduction 

In this section we describe the charged Higgs boson discovery potential of the ATLAS detec- 
tor in the (m//±, tan (3) parameter space which has been investigated using the ATLFAST 



and PYTHIA 5.7 p4| simulation packages. This is a particle-level simulation performed at 
^/s = 14 TeV, but with the detector resolutions and efficiencies parametrized from full de- 
tector simulations. It is assumed that the mass scale of supersymmetric partners of ordinary 
matter is above the charged Higgs bosons so that decays into supersymmetric partners 
are forbidden A central value 175 GeV is used for the top-quark mass. 



The decays if^ tb and tu are the dominant channels in most of the parameter 

space |T^. The decay channel tu has been studied extensively for ATLAS for 

mjyi < rrit, and the signal appears as an excess of r leptons [^. The channel — »• Wh^ 



is only relevant in a tiny range of MSSM parameter space but it constitutes a unique test for 
MSSM and may be sensitive to the singlet extension to MSSM, i.e., NMSSM. The cs 
channel is studied as a complement to the r-lepton channel: if the charged Higgs is detected 
by observing the excess of r-leptons over the SM prediction, then the cs channel could be 
used to measure mH±. Discovery is possible through the — > tb channel for low ( < 3) and 
large ( > 25) tan /? values up to masses mH± ~ 400 GeV. In the following, a brief description 



of the analysis is presented; details can be found elsewhere [32 



5.2 Discovery Potential 

(i) t — s> bH^ — s> bcs, mH± < mt: ti events are generated through gg, qq ti with one 
top-quark decaying into the charged Higgs, and the other into W, t — >■ Wb Ivb. The major 
background is tt production itself with both top-quark decaying into VT's; one of the VT's 
goes to jets and the other to leptons. This process is studied for tan/? = 1.5 and = 110 
and 130 GeV. The events with a final state consisting of two b-tagged jets (|r/| < 2.5, and 
Pt > 15 GeV), and a single isolated lepton (|?7| < 2.5, > 20 and pj. > 6 GeV) are 
selected and the charged mass peak is searched for the di-jet mass distribution rrijj. The 
combinatorial background is reduced by applying a b-jet veto and a jet-veto on extra jets. 
Fig. 1^ shows the di-jet mass distribution for both the signal and the background. This 
channel complements the — > tu channel in that if the is detected by observing the 
excess of r-leptons, the cs channel can be used to determine mH±. 

(ii) t — > bH''', W*h°, mH± < nit: The production mechanism is the same as in 

the previous case, but here, — * W*h^ , with hP — > bb. The final state contains two VF's, 
one of which is off-shell and one of which decays to leptons and the other to jets. The major 
backgrounds are tibb and tiqq followed by the decays of the top-quarks as described above. 
The present channel is studied for mH± = 152 GeV and for tan/3 = 2 and 3 corresponding 
to m^o = 83.5 and 93.1 GeV respectively. We search for an isolated lepton, four b-tagged 
jets (py > 30 GeV) and at least two non b-jets with pip > 30 GeV. The details of this 



analysis can be found in It suffices to say that although the backgrounds are over 
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Figure 9: For the cs channel, the expected rrijj distribution from signal and background 

events (solid) and from the background (dashed) for mff± = 130 GeV and tan/3 = 1.5 and for an 
integrated luminosity of 30 fb~^. Errors are statistical only. 

two orders of magnitude higher that the signal at the start, we propose a reconstruction 
method which permits the extraction of the signal with a significance exceeding 5a in the 
low tan/? (1.5 — 2.5) region. At high tan/5, though the reconstruction remains comparable 
the signal rate decreases so significantly that discovery potential vanishes in this region. 
Fig. |10| shows the charged Higgs mass reconstruction for tan j3 = 2. 

(iii) mH± > nit: Above the top-quark mass, we consider the production of through 
the 2 — * 2 process gb tH^. Two decay channels of are examined in details, H"^ th 
and —>■ Wh^ —>■ Wbb. In both cases the major background comes from ttb and ttq 
events. In either case, we search for an isolated lepton, three b-tagged jets and at least 
two non b-jets. The details of these analyses can be found elsewhere p2| , |33| . Discovery 
is possible through the tb channel for low (< 3) and for high (> 25) tan/5 up to 

400 GeV H. Fig. |ll] shows the charged Higgs mass reconstruction for tan/5 = 1.5 
300 GeV. On the other hand, the — >■ Wh'^ channel presents no discovery 



and mH± 

potential for the charged Higgs in the MSSM. Initially, the total background is at least three 
orders of magnitude higher than the signal in the most favorable case studied (tan/? = 3). 
We propose a reconstruction technique which improves the signal-to-background ratios by 
two orders of magnitude. However, this improvement is still not enough to observe a clear 
signal; for example, at tan/5 = 3, a significance of only 3.3 can be expected after three years 



of high luminosity operation 



5.3 Conclusions 

The possibility of detecting the charged Higgs through the decay channels cs, 

Wh^, and tb with the ATLAS detector has been studied as a function of tan /5, below 
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Figure 10: For the — > W*h^ channel, the reconstructed mass distribution from 
signal+background events (sohd) and from background events (dashed) for = 152 GeV, 

tan/3 = 2, and for and integrated luminosity of 300 fb~^. Errors are statistical. 



> 

o 

(N 

£h 

D 

B 
x> 

o 

I 

o 
o 

o 

> 
o 



> 



150 



100 



50 — 



200 




Figure 11: Signal and background distributions for the reconstructed invariant mass m^f, for mjj± 
300 GeV, tan/3 = 1.5 and an integrated luminosity of 30 fb^^. Errors are statistical only. 



and above the top-quark mass. Below the top-quark mass and at low tan (3, both channels 
— >■ cs and Wh^ present significant discovery potential. These two channels would 

complement the tu searches in that if the latter is observed through the excess of 

r-leptons, the former channels can be used to measure the mass of the charged Higgs. Above 
the top-quark mass, the process th presents a significant discovery potential in the 

low and the high tan (3 regions up to 400 GeV. 
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Light stop effects and Higgs boson searches at the LHC. 

G. Belanger, F. Boudjema, A. Djouadi, V. Ilyin, 
J.L. Kneur, S. Moretti, E. Richter-Was and K. Sridhar 

Abstract 

We analyze the effects of light top squarks with large mixing on the search of the 
lightest Higgs boson of the Minimal Supersymmetric extension of the Standard Model 
at the LHC. We discuss both the stop loop effects in the main production and decay 
processes, and the associated production of top squarks with the lightest Higgs boson. 

1 Introduction 

The third generation fermions, and especially the top quark because of its large Yukawa 
coupling, play an important role in the mechanism of electroweak symmetry breaking and 
the properties of the Higgs bosons Recall that if the top quark were rather light, the 
Minimal Supersymmetric extension of the Standard Model (MSSM) would have been already 
discarded since the lightest Higgs boson h that it predicts would have been lighter than the 
Z boson, Mh < Mz and would have not escaped detection at LEP2. The contribution of 
the top quark and its SUSY partners to the radiative corrections to can push the mass 
value up to Mh ~ 135 GeV ||2|, beyond the reach of LEP2. The mixing in the stop sector 
is also important since large values of the mixing parameter At = At + ii/ tan/3 [where At 
is the trilinear coupling, /i the higgsino mass parameter and tan/5 the ratio of the vev's of 
the two Higgs doublets which break the electroweak symmetry; see Ref. for the SUSY 
parameters] can increase the h boson mass for a given value of tan/3 [0. 

On the other hand, while the sfermions of the two first generations can be very heavy, 
naturalness arguments suggest that the SUSY particles that couple substantially to the Higgs 
bosons [stops, sbottoms for large tan/5, and the electroweak gauginos and higgsinos] could 
be relatively light. In this respect, the case of the stop sector is special: because of the large 
nit value, the mixing in this sector can be very strong, leading to a mass eigenstate ti lighter 
that all other squarks, and possibly lighter than the top quark itself. At the same time, 
again because of the large mixing, this particle can couple very strongly to the MSSM Higgs 
bosons and in particular to the lightest CP-even particle h. 

At the LHC, a light stop with large couplings to Higgs bosons can contribute to both 
the h production in the main channel, the gluon-gluon fusion mechanism gg — > h, and to 
the main detection channel, the two-photon decay h —>■ 77. The effects can be extremely 
large, making this discovery channel possibly useless at the LHC [4-6]. On the other hand, 
because of the enhanced couplings and phase-space, associated production of stops and the 
h boson at the LHC, pp qq/gg — ^ iiiih, might have sizeable cross sections [7-10]. 

It is thus crucial to investigate how and when this scenario occurs and what other conse- 
quences then follow at the LHC. The purpose of our working group contribution is to update 
and complement the various analyses [5-10] which have been made on this subject. 
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2 Stop parameters and phenomenological constraints 



We start our discussion by recalling the parameters that define the stop masses, mixing angle 
and the titih coupling. The stop mass eigenstates are defined through the mixing angle 6];, 
with the lightest stop ti being ti = cos 6]; ti — smOi t^. With the effective trilinear mixing 
parameter, At = At + ^/ tan/3, one has for the masses and the mixing angleQ 

tan(2%) = -TT. YTT^ l^TTT. s 9 n sin(26'i) = — 5 5- (1) 

^ ^' m\ -m% + iMlcos2/?(l - fsL) ^ *^ m? - m? ^ ^ 



*1,2 



m. 



+ 2 



(2) 



where '^Qg, cire the soft-SUSY breaking scalar masses and the dots stand for the D 

terms oc M| cos 2/3. Note that in order to enhance the mixing, sin(2^^j) ~ 1, one needs to 
make At large and/or have the soft-SUSY masses almost equal: fn 



The titih 



vertex writes 



iiiih 



+ 



-9 



rrit cos a 



Mw sin (3 



(y4t — yutana) sin6'j cos^ — rrit 



M| sin (3 . 



rrit cos a 



sin(a + /?) 



9 



Ml 



w 



12 2 
[2^3 si^^ ^w) cos^ Oi + - sin^ 9w sin^ 9i 



\ 2 

rriT ) + rrij 



(3) 



where in the last line we neglected the D-term contributions and assumed the limit of large 
Ma to be in the decoupling regime. As can be seen, in the presence of large mixing with 
large splitting between the two stop eigenstates, the iiiih coupling can be particularly large. 
In the case of no mixing, only the top contribution survives and the coupling titih is of the 
order of tih coupling. Taking this limit as a reference point, the strength of the titih vertex 

r n 2 

can be normalized through Rt^ = MwVtj^fj^/ (gm^) . 

We now summarize the constraints which can be imposed on the stop parameters: 

• The model independent mass limit on the lightest stop is obtained from direct searches 
at LEP, m^-^ > 90 GeV |TT|]. However, if the ti and the xo LSP are not too close in 
mass, a stronger limit, m^^ > 120 GeV [Q, is available from Tevatron analyses. For 
bottom squarks, a limit > 250 GeV is available from Tevatron data in the case of 
no-mixing |T^. 

• If stops are too light, the radiative corrections to the h boson mass are not large enough 



and the limit > 90 GeV from LEP searches plays an important role. 



^The sign conventions for At here is opposite to the one adopted in Refs. ||] and [Q. Accordingly, the 
sign convention for the mixing angle is opposite to the one of Ref. M where ii — cosOj: <l + sinS^ iji. 
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As in the case of top/bottom splitting in the Standard Model, the stop/sbottom dou- 
blet can contribute significantly to electroweak precision observables through the p 
parameter. In particular, if stops strongly mix and have large couplings, the contribu- 
tions to Ap can exceed the value Ap < 0.0013 imposed by data 



Some values of the stop parameters might induce color and charge breaking minima 
(CCB). Since the naive constraints based on the global minima may be too restrictive, 
we will take into account the tunneling rate [for wide range of parameters, the global 
CCB minimum becomes irrelevant on the ground that the time required to reach the 
lowest energy state exceeds the present age of the universe], which leads to a milder 



constraint which may be approximated by [|T^]: + < 7.5(M?^ + M?. ] 




Figure 1: (a) Constraint from Ap < 0.0013 (full line), Mj, > 90 GeV (dash-dot), CCB (dash) and 
(dash) for tan/3 = 10, /i = 400 GeV, m^-^ = 120 GeV and Ma = I TeV; the Mh constraint 
for tan/3 = 2.5 is also shown (dot), (b) Equipotential lines (dotted) for the normalized coupling 
Rl_^ = 1,10,50,100 with tan/3 = 10 and /i = 400 GeV. The exclusion regions corresponding to 
Ap < .0013 and < 90 GeV are also reproduced. 

Fig. 1 shows how the parameter space is restricted by the previous constraints and which 
values of the ratio Ri^ are allowed. In Fig. la, the excluded region in the plane (cos 6'^ , mi^ ) is 
within the respective boundaries indicated. Note that for cos6'j ~ 1, the Ap constraint also 
excludes the region to the right of the second branch of the Ap curve where the present limit 
on the mass of the sbottom is contained. Requiring m^^ > 250 GeV excludes the region to 
the right of the curve. The CCB constraint for p = 800 GeV is also displayed, the excluded 
region lies between the two "CCB, p = 800" curves. In Fig. lb, we show the equipotential 
lines for the normalized coupling Ri^. The exclusion regions corresponding to Ap < .0013 
and Mh < 90 GeV are also reproduced. In all cases M2 = —fj. and a common gaugino mass 
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at the GUT scale are assumed. Note that one has to make sure that the hghtest stop is not 
the LSP, as has been always verified in our analysis. Considering that the CCB constraint 
is rather uncertain, it is also worth pointing out that the one used in our analysis hardly 
precludes points which are not already excluded by the Ap and mh constraints. 



3 Higgs boson signals at the LHC 

In this section, we will discuss what might happen to the search for the lightest MSSM Higgs 
boson h at the LHC, if one allows all sparticles but the stops (and to a lesser extent the 
charginos and neutralinos) to be rather heavy. We will first discuss the effects of stop loops 
in the gluon-gluon fusion mechanism, gg h, and in the main Higgs detection channel, the 
two-photon decay h — > 77, and then discuss the associated production of stops with the 
light Higgs boson h and possibly A. 



3.1 Stop loop effects 



2 



Since the htiti vertex eq. (2) does not have a definite sign [for no mixing the positive 
component dominates while for maximal mixing the negative component | sin^(2^£)(m?^ — 
m?^) is the leading one], the stop loop contributions can interfere either destructively or 
constructively with the top loop contributions in the gg —>■ h and h 'y-f processes. Noting 
that while for gg h only top/stop loops are present, for the decay h — > 77, the additional 
contributions from W loops are dominant and have a destructive interference with the top 
contributions. This means that if the rate for h ^ gg is suppressed, there will be a slight 
increase in /i ^ 77 decay width and vice versa. Therefore either the rate for the inclusive 
channel gg —>■ h 'y-f is enhanced or the rate for the associated Higgs production pp 



Wh, Zh,tih IjT^ with — > 77 is enhanced. It is important to stress that, in any case, the 
rate for the associated tih production with the subsequent decay h —>■ bb is hardly affected 
by stop loops and will always help in these scenarii, as will be discussed later. 

We begin our analysis by defining the ratio R^^ = which is the branching ratio of 

the lightest SUSY Higgs boson decay into two photons over that of the SM for the same Higgs 
mass. In the decoupling regime. Ma ^ Mz, this ratio is affected only by SUSY-particle 
loops; in this case the ratio is also sensibly the same as the ratio for associated production 
of the h boson with W, Z bosons and/or with tt pairs, with h decaying into 77. We also 
define Rgg-yy as the ratio for the signal in the direct production channel gg h times the 
branching ratio for the h ^ y'j decay in the two models. The gg and 77 decay widths are 
obtained^ with the help of the program HDECAY ||16|| . 

Fig. 1^ summarizes the contribution of stop loops to these ratios, for tan/5 = 2.5, /i = 
— M2 = 250 GeV and Ma = 1 TeV. To maximize the effect of stop mixing, sin(26'j) ~ 1, we 
assume that rrig^ ~ rhjy^^. From this figure, one can see that: 

^Note that the ratios of gg decay widths and production cross sections are almost the same: large QCD 
corrections cancel out in the ratios when the dominant contribution comes from the top loops, and the 
corrections to the top and stop contributions are practically the same; see Ref. 
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Figure 2: Higgs boson (/i) production and decay ratios at the LHC for rFig^ ~ fhfj ^ at tan /3 = 2.5 
and large Ma- Figures are scanned over and At within the constraints discussed above. 
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- The /i — s> 77 branching ratio is only mildly affected [less than ~ 30%] by the contribu- 
tions of the stop loops which can be of either sign. This is mainly due to the fact that the 
W contribution to the /177 vertex is largely dominant in the decoupling limit. 

- The hgg coupling is always reduced compared to the SM case for large stop mixing 
and rather light stops can lead to a strong reduction in the rate of the inclusive production 
channel gg h. The suppression factor can be as low as 1/10 whereas a benchmark for 
discovery is about ~ 1/2 [although this benchmark depends slightly on the Higgs boson 
mass]. The suppression occurs for rather large, though not maximal, Higgs boson masses 
where the efficiencies are better than for smaller Higgs masses. 

- For very heavy stops which should decouple from the hgg and /177 vertices, the ratio 
Rgg-y-y couM bc diffcreut from unity since charginos could be also hght and might give small 
contributions to the /i ^ 77 decay width in the MSSM. 



3.2 Associated Higgs production with stops 

If the mixing in the stop sector is large, one of the top squarks can be rather light and at 
the same time, its couplings to the Higgs boson can be strongly enhanced. The associated 
production process pp — > qq/gg —>■ htiti might then be favored by phase space and the cross 
sections might be significantly large. This process is thus worth investigating at the LHC 

In view of the implementation of the process pp qq/gg titih into an event generator, 
it is useful to give a "model independent" description of the production cross section in the 
continuum, in terms of the parameters m^^, [besides mt etc...]. One can tabulate, in a 
way which can be read externally, the cross section according to selected valuesQof Mh,m^_^ 
together with the coupling V^j^fi [for simplicity and as a first step, one can take the vertex 
Vi^l^fi such that i?^^ = 1, i.e. in the large M4 limit, no t mixing and D-terms]. 

The generator of partonic events for pp titih can be created by using the package 
CompHEP and may be down-loaded at this http address [^. The events can be used as 



an external process input in PYTHIA [2C] or ISA JET for further decay and hadronization to 
simulate full events at the level of detectable particles. The titih coupling is evaluated as 
a user's function thus allowing for an interface with any SUSY model. The generator also 
includes, as an option, the event generation of the SM process pp — > qq/gg tt+Higgs. 

As an illustration, defined reference cross section^ calculated with the help of CompHEP 
are displayed in Fig. ^ The cross sections are shown as functions of M/i(mjJ for given 
values of mi^{Mh), for a titih vertex in the limit of large Ma, no mixing and no D-terms, 
as discussed above. Also shown are the cross sections for the processes pp — >■ tth, ttZ [where 
only the dominating contributions of the gg initiated subprocesses are included] and iiiiZ 
[where the vertex has been computed with cos^ 6^ = 1/2, i.e. maximal mixing, and has to be 
rescaled by a factor (cos^ 9i/2 — 2/3s^) for other mixing values,]. We have used the CTEQ4 
structure functions |p2| with a scale set at the invariant mass of the subprocess. 



^Of course, in reality, the situation is slightly more complicated since the two masses mi_^,Mh and the 
coupling Vt-^^tih depend on the mixing and are thus inter- related 
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Note that the complete analytical expressions of the pp gg/qq OT+Higgs are given in Ref. [0. 
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Figure 3: The cross section pp titih (and similar processes) at the LHC as functions of m^^ 
(left) and (right) for a range of and m^^ values. See text for details. 

As can be seen, the pp — > titih cross section can be large for small values of the stop 
and the Higgs masses, but drops precipitously with m^^ and to a lesser extent with Mh- The 
cross section is more than order of magnitude larger than a{pp titiZ) and can exceed the 
one for the SM-like process pp —>■ tih for strong enough mixing Ri ^ 1 and light ti. 

If one takes the value a{titih) > 300 fb as a benchmark cross section value for observing 
this process at the LHC, and using the constraint on the maximum values of the titih 
coupling, values of m^-^ > 250 GeV are hardly accessible at the LHC. This is shown in Fig. ^ 
where the pp titih cross section is shown as a function of m^^ taking all soft squark masses 
equal for tan/3 = 2.5 and imposing Ap < .0013. A scan on the common soft breaking scalar 
mass and At has also been performed; shown are points that pass the criteria (Jij^h > 300 fb 
and for which m^^ > 150 GeV. Larger values of the stop mass can be reached if ones allows 
a 2a variation on the Ap constraint, as shown in the figures at the bottom. 



3.3 Comparison of inclusive and associated production processes 

Let us now make a global discussion on the stop effects in both type of processes for Higgs 
boson production, gg ^ h —>■ j'-f and pp titih. The two cross sections are shown in 
Fig. ^ in the decoupling limit for tan/5 = 2.5 and equal soft breaking scalar masses. As can 
be seen, the suppression of the rate in the inclusive production channel is compensated by 
a rate increase in the associated production channel. When the suppression factor in the 
inclusive production is below 0.5 making discovery in this channel difficult, the cross section 
for the process titih is above 200 fb. As discussed previously, a benchmark value for the 
cross section allowing the discovery of the Higgs boson in the pp titih channel has been 
estimated to cr ~ 300 fb. Therefore for some values of the parameters, neither the inclusive 
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Figure 4: The pp titih cross sections at the LHC as a function of m^^ (left) and a scan on itIq 
and At (right). Shown also are points that pass cr^^f^/i > 300 fb and m^-^ > 150 GeV, imposing 
Ap < .0013 (top) or Ap < .0026 (bottom). 
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nor the pp titih channels can be accessed. However, one also sees that for these same 
points one can without difficulty use the usual pp Wh/Zh,tih search modes. 
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Figure 5: The cross sections for the inclusive and associated Higgs production at the LHC for 
~ "^c/3ji tan/3 = 2.5 and large Ma- Figures are scanned over and Af. 

Therefore with the remark that the process pp — » tih [with h bb] should allow for 
Higgs boson discovery at the LHC within this scenario, one should salvage the detection the 
h boson with the bonus that the stop should also be observed. Even though one may have to 
wait for the higher luminosity stage, the scenario with light stops and large couplings offer 
much better prospects than previously thought. 

The assumption of an equal value for the soft scalar masses at the weak scale is rather 
unnatural [see later in mSUGRA] and could be relaxed. To illustrate the fact that large 
suppression factors in the inclusive production channel, though not as dramatic as in the 
previous case, stiU occur we show in Fig. ^ typical R ratios for unequal values of the soft 
masses. What is most interesting is that, as soon as sm{26^) ^ 1, the non-diagonal decay 
channel ti iih opens up and can have an appreciable branching ratio. This can be seen 
by inspection of the Vi^i^^^ coupling, for which the leading component is proportional to: 
Hii2/i ^ qK^^w) sin 46*^(771?^ — m|). Considering that if the ^2 mass is not excessively large, 
^2 is produced in abundance and this cascade decay can provide more Higgs bosons than 
through the continuum pp iiiih production. 

Perhaps even more interesting, is the case when is not too large. For large values 
of At, and even when sin2^f ~ 1, one can have a large decay rate ^2 ti^ since the 
Atit2 coupling can be large H^JjA o*^ g^t/ {'^Mw){At/ tan/5 — /i), as shown in Fig. |^. This 
coupling is generally larger than the t2tiH coupling and hence, within these scenarii, the 
decay t2 —>■ tiA is most likely to occur than the decay into the heavier H boson, t2 — »■ Hti. 



Finally, let us make a few comments on the case of the minimal Supergravity model [^ , 
where the only input parameters are the universal scalar mass mo, the universal gaugino 
mass parameter mi/2, the trilinear coupling Aq, tan/? and the sign of the /i parameter. The 
parameters mo,mi/2 and Aq are chosen at the GUT scale and their evolution down to the 
weak scale is given by the RGB's [03. Proper breaking of the electroweak symmetry is also 



assumed, which fixes the parameter In what follows, the RGB's and the proper BW 



symmetry breaking are solved using the program SUSPECT ^5 . 
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Figure 6: Higgs boson production rates at the LHC for unequal soft breaking scalar masses in 
decoupling limit Ma =lTeV and a{pp — > i2iih) vs cr{pp — > iiiih). 
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Figure 7: Higgs boson production and decay rates at the LHC for low values of Ma- 
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In the mSUGRA case the cross section can be as large as in the case of the unconstrained 
MSSM, but in a relatively smaller area of the SUSY parameter space. This is essentially 
due to the fact that it is generically very difficult to have almost degenerate II and in 



mSUGRA, so that the stop mixing angle which is controlled by the ratio At/{'mi — 



can become large only for very large At. Moreover in the RG evolution l^^l tends to 
decrease when the energy scale is decreasing from GUT to low-energy. This makes a large 
At value at low energy less likely, since Aq = ^^(GUT) would have to be even larger, which 
may conflict with e.g. the CCB constraints. The only way to have an increasing \At\ when 
running down to low energy is if < with Aq small enough. This requires a large mi/2 
value, which implies not too small m^-^ . 

Thus the mixing in the stop sector is, in general, not as large as in the unconstrained 
MSSM and the titih coupling for instance is, in general, smaller than in the previous case. 
This implies that the rate for the inclusive production and detection channel gg h ^ '-fj 
[in the decoupling limit] is not as dramatically different from the rate in the SM, as it can 
be in the unconstrained MSSM. Furthermore, the milder mixing results in a smaller cross 
sections for the process pp —>■ titih as is shown in Fig. 8 for LHC energies. However, for 
large tan (3 values, the pseudoscalar A boson tends to be rather light in mSUGRA, opening 
the possibility for the decay t2 tiA to occur with an appreciable rate as shown in Fig. 8b. 

Note that one should also easily observe the pseudoscalar A boson in the loop mediated 
process gg ^ A since the rate is strongly enhanced for large tan/3 and, because of CP- 
invariance, light stop [or sbottom] loops cannot contribute to the process. 
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Figure 8: Cross sections in mSUGRA at LHC: a) a{pp ->■ iiiih) for m^/a = 0.3 TeV, Aq = 2 TeV, 
tan/? = 2.5,30. b) a{pp iihhJihhJihA) for mo = 0.2 TeV, = 0.3 TeV and tan/? = 35. 
For the spectrum, SUSPECT is used in a) and ISAJET in b). 

References 



64 



[1] For a review of the Higgs sector in the MSSM, see J.F. Gunion, H.E. Haber, G.L. Kane, 
S. Dawson, The Higgs Hunter's Guide, Addison-Wesley, Reading 1990. 

[2] For a recent analysis, see M. Carena et al., hep-ph/0001002. 

[3] For a review on SUSY see e.g.: H.E. Haber, G.L. Kane, Phys. Rep. 117, 75 (1985). 

[4] B. Kileng, Z. Phys. C63 (1993) 87; B. Kilcng, P. Osland, P.N. Pandita, Z. Phys. C71 
(1996)87 and hcp-ph/9506455; G. L Kane, G. D. Kribs, S.P Martin and J. D. Wells, 
Phys. Rev. D50 (1996) 213. 

[5] A. Djouadi, Phys. Lett. B435 (1998) 101; ibid. Mod. Phys. Lett. A14 (1999) 359. 

[6] G. Belanger, F. Boudjema and K. Sridhar, hep-ph/9904348. 

[7] A. Djouadi and J. L. Kneur and G. Moultaka, Phys. Rev. Lett. 80 (1998) 1830; zbzd 
hcp-ph/9903218. 

[8] G. Belanger, F. Boudjema, T. Kon and V. Lafage, hep-ph/9711334. 

[9] A. Dedes and S. Moretti, Phys. Rev. D60 (1999) 015007; ibidEur. Phys. J. CIO (1999) 
515; ibid hep-ph/9909526. 

[10] G. Belanger, F. Boudjema, T. Kon and V. Lafage, hep-ph/9907207. 

[11] For an update on the limits on the SUSY Higgs bosons and particle masses, see the 
talks of the LEP collaborations for the LEPC Committee, Nov. 1999: 
A. Blondel, ALEPH talk: 

http://alephwww.cern. ch/ALPUB/seminar/lepc_nov99/lepc.pdf ; 

J. Marco, DELPHI talk 

http://delphiwww.cern.ch/ of f line/physics_links/lepc.html. 
G. Rahal-Callot, L3 talk 

http : //13www . cam . ch/conf erences/ps/RahalCallot_LEPC_L3_Nov99 . ps . gz. 

P. Ward, OPAL talk 

http : // wwwl . cern . ch/0pal/talks/pward_lepc99 . ps . gz. 

[12] T. Affolder et al, The CDF Collaboration, FERMILAB-PUB-99/311-E. 

[13] See for instance, G. Altarelli, hep-ph/96 11239; J. Erler and P. Langacker, hep- 
ph/9809352; G.C. Cho et al, hep-ph/9901351; A. Djouadi, hep-ph/9911468. 

[14] A. Kusenko, P. Langacker and G. Scgre, Phys. Rev. D54 (1996) 5824; A. Kusenko and 
P. Langacker, Phys. Lett. 391 (1997) 29; For a summary of "naive CCB" constraints, 
see J. A. Casas, hep-ph/9707475. 

[15] A. Djouadi, M. Spira, P.M. Zerwas, Phys. Lett. B264 (1991) 440; S. Dawson, Nucl. 
Phys. B359 (1991) 283; M. Spira et al, Nucl. Phys. B453 (1995) 17; S. Dawson, A. 
Djouadi, M. Spira, Phys. Rev. Lett. 77 (1996) 16. 



65 



[16] A. Djouadi, J. Kalinowski and M. Spira, Comp. Phys. Comm. 108 (1998) 56. 

[17] For an update of these cross sections at the LHC, see: Z. Kunszt, S. Moretti, W. Stirhng, 
Z. Phys. C74 (1997) 479; M. Spira, Fortsch. Phys. 46 (1998) 203. 

[18] For a description of CompHep, see A.Pukhov et al, "CompHEP user's manual, v3.3". 
Preprint INP MSU98-41/542, 1998; hep-ph/9908288 

[19] /af s/cern. ch/cms/physics/PEVLIB/susy-tth. 

[20] T. Sjostrand, PYTHIA 5.7 and JETSET 7.4 Physics and Manual, hep-ph/9508391. 
[21] ISAJET 7.48, H. Baer, F. E. Paige, S.D. Protopescu and X. Tata, hep-ph/0001080. 
[22] CTEQ Collaboration, Phys. Rev. D51 (1995) 4763. 

[23] A.H. Chamseddine, R. Arnowitt and P. Nath, Phys. Rev. Lett. 49 (1982) 970; R. 
Barbicri, S. Fcrrara and C.A Savoy, Phys. Lett. B119 (1982) 343; L. Hall, J. Lykken 
and S. Weinberg, Phys. Rev. D27 (1983) 2359. 

[24] See for instance, D.J. Castaiio, E.J. Piard and P. Ramond, Phys. Rev. D49 (1994) 
4882; W. de Boer, R. Ehret and D.I. Kazakov, Z. Phys. C67 (1994) 647; V. Barger, 
M.S. Berger and P. Ohmann, Phys. Rev. D49 (1994) 4908; M. Drees and S. Martin, 
hep-ph/9504324. 

[25] A. Djouadi, J.L. Kneur and G. Moultaka, hep-ph/9901246; the program can be found 
at: http://www.lpin. univ-montp2.fr: 7082/ djouadi/GDR/mssm4.htiiil. 



66 



Double Higgs production at TeV Colliders 
in the Minimal Supersymmetric Standard Model 

R. Lafaye, D.J. Miller, M. Muhlleitner and S. Moretti 

Abstract 

The reconstruction of the Higgs potential in the Minimal Supersymmetric Standard 
Model (MSSM) requires the measurement of the trilinear Higgs self-couplings. The 
'double Higgs production' subgroup has been investigating the possibility of detecting 
signatures of processes carrying a dependence on these vertices at the Large Hadron 
Collider (LHC) and future Linear Colliders (LCs). As reference reactions, we have 
chosen gg — > hh and e~^e~ — > hhZ, respectively, where h is the lightest of the MSSM 
Higgs bosons. In both cases, the Hhh interaction is involved. For niH ^ 2m/i, the 
two reactions are resonant in the H ^ hh mode, providing cross sections which are 
detectable at both accelerators and strongly sensitive to the strength of the trilinear 
coupling involved. We explore this mass regime of the MSSM in the /i — > 66 decay 
channel, also accounting for irreducible background effects. 

1 Introduction 

Considerable attention has been devoted to double Higgs boson production at future e^e~ 
and hadron colliders, both in the Standard Model (SM) and the MSSM [0, |, I). For the 
SM, detailed signal-to-background studies already exist for a LC environment 0, for both 
'reducible' and 'irreducible' backgrounds 0, |[, which have assessed the feasibility of experi- 
mental analyses. At the LHC, since here the typical SM signal cross sections are of the order 
of 10 fb [0], high integrated luminosities would be needed to generate a statistically large 
enough sample of double Higgs events. These would be further obscured by an overwhelm- 
ing background, making their selection and analysis in a hadronic environment extremely 
difficult. Thus, in this contribution we will concentrate only on the case of the MSSM. 

In the Supersymmetric (SUSY) scenario, the phenomenological potential of these reac- 
tions is two-fold. Firstly, in some specific cases, they can furnish new discovery channels for 
Higgs bosons. Secondly, they are all dependent upon several triple Higgs self-couplings of 
the theory, which can then be tested by comparing theoretical predictions with experimental 
measurements. This is the first step in the reconstruction of the Higgs potential itsellQ. 

The Higgs Working Group (WG) has focused much of its attention in assessing the 
viability of these reactions at future TeV colliders. However, the number of such processes is 
very large both at the LHC and a LC P, |^, so only a few 'reference' reactions could be studied 
in the context of this Workshop. Work is in progress for the longer term, which aims to cover 
most of the double Higgs production and decay phenomenology at both accelerators 

^^The determination of the quartic self- interactions is also required, but appears out of reach for some 
time: see Refs. ^, |^ for some cross sections of triple Higgs production. 
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These reference reactions were chosen to be the gluon-fusion mechanism, gg hh, for 
the LHC (see top of Fig. |l]) and the Higgs-strahlung process, e+e~ hhZ, for the LC (see 
bottom of Fig. |^), where h is the hghtest of the MSSM scalar Higgs bosons. The reason 
for this preference is simple. Firstly, a stable upper limit exists on the value of rrih, of the 
order of 130 GeV, now at two-loop level [0, so that its detection is potentially well within 
the reach of both the LHC and a LC. In contrast, the mass of all other Higgs bosons of 
the MSSM may vary from the electroweak (EW) scale, 0{mz), up to the TeV region. In 
addition, as noted in Ref . 0] , the multi-fe final state in gg ^ hh ^ bbbb, with two resonances 
and large transverse momenta, may be exploited in the search for the h scalar in the large 
tan f3 and moderate region. This is a corner of the MSSM parameter space that has so far 
eluded the scope of the standard Higgs production and decay modes p. (The symbol A here 
denotes the pseudoscalar Higgs boson of the MSSM, and we reserve the notation H for the 
heaviest scalar Higgs state of the model.) However, this paper will not investigate the LHC 
discovery potential in this mode, given the very sophisticated treatment of the background 
(well beyond the scope of this note) required by the assumption that no h scalar state has 
been previously discovered (see below). This will be done in Ref. Furthermore, the 
gg hh and e+e^ hhZ modes largely dominate double Higgs production 0, ^, at least 
for centre-of-mass (CM) energies of 14 TeV at the LHC and 500 GeV in the case of a LC, 
the default values of our analysis. (Notice that we assume no polarization of the incoming 
beams in e~^e~ scattering.) Finally, when mn ^ 2m/i, the two reactions are resonant, as 
they can both proceed via intermediate states involving H scalars, through gg ^ H and 
e'^e^ HZ, which in turn decay via H ^ hh Thus, the production cross sections 
are largely enhanced 0] (up to two orders of magnitude above typical SM rates at the 
LHC 0) and become clearly visible. This allows the possibility of probing the trilinear Hhh 
vertex at one or both these colliders. 

The dominant decay rate of the MSSM h scalar is into bb pairs, regardless of the value 
of tan/5 P|. Therefore, the final signatures of our reference reactions always involve four 
6-quarks in the final state. (In the case of a LC environment, a further trigger on the 
accompanying Z boson can be exploited.) 

If one assumes very efficient tagging and high-purity sampling of 6-quarks, the background 
to hh events at the LHC is dominated by the irreducible QCD modes ||10|- Among these, 
we focus here on the cases qq, gg bbbb, as representative of ideal 6-tagging performances. 
These modes consist of a purely QCD contribution of O^a^), an entirely EW process of 
C(ag^) (with no double Higgs intermediate states) and an 0{a'^al^) component consisting 
of EW and QCD interactions. (Note that in the EW case only qq initiated subprocesses are 
allowed at tree-level.) For a LC, the final state of the signal is bbbbZ, with the Z reconstructed 
from its decay products in some channel. Here, the EW background is of 0{al.^) away from 
resonances (and, again, contains no more than one intermediate Higgs boson), whereas the 
EW/QCD background is proportional to (afog^). 

In general, EW backgrounds can be problematic due to the presence of Z vectors and 
single Higgs scalars yielding bb pairs, with the partons being typically at large transverse 
momenta and well separated. In contrast, the QCD backgrounds involve no heavy objects 
decaying to bb pairs and are dominated by the typical infrared (i.e., soft and collinear) QCD 
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gg to double Higgs fusion at the LHC: gg hh 
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double Higgs-strahlung at a LC: e'^e hhZ 
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Figure 1: Diagrams contributing to gg hh (top) and e^e hhZ (bottom) in the 
MSSM. 

behavior of the partons in the final state. However, they can yield large production rates 
because of the strong couplings. 

In this study, we investigate the interplay between the signal and background at both 
colliders, adopting detector as well as dedicated selection cuts. We carry out our analysis 
at both parton and hadron level. The plan of this note is as follows. The next Section 
details the procedure adopted in the numerical computation. Sect. displays our results 
and contains our discussion. Finally, in the last section, we summarize our findings and 
consider possible future studies. 



2 Calculation 

For the parton level simulation, the double Higgs production process at the LHC, via gg fu- 
sion, has been simulated using the program of Ref. to generate the interaction gg hh, 



with the matrix elements (MEs) taken at leading-order (LO) for consistency with our treat- 
ment of the background. We then perform the two h bb decays to obtain the actual 
46-final state. For double Higgs production at a LC, we use a source code for the signal 
derived from that already used in Ref. 0. At both colliders, amplitudes for background 
events were generated by means of MadGraph |jl2| and the HELAS package Note that 



interferences between signal and backgrounds, and between the various background contri- 
butions themselves, have been neglected. This is a good approximation for the interferences 
involving the signal because of the very narrow width of the MSSM lightest Higgs boson. 
Similarly, the various background subprocesses have very different topologies, and one would 
expect their interferences to be small in general. 

The Higgs boson masses and couplings of the MSSM can be expressed at tree-level 
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in terms of the mass of the pseudoscalar Higgs state, m^, and the ratio of the vacuum 
expectation values of the two neutral fields in the two iso-doublets, tan (3. At higher order 
however, top and stop loop-effects can become significant. Radiative corrections in the one- 
loop leading approximation are parameterized by 



V2n^ sin^ 



rric 



13 



log^ 



where the SUSY breaking scale is given by the common squark mass, ms, set equal to 1 TeV 
in the numerical analysis. If stop mixing effects are modest at the SUSY scale, they can 
be accounted for by shifting m| in e by the amount Am| = 74^[1 — A^/(12m|)] {A is the 
trilinear common coupling). The charged and neutral CP-even Higgs boson masses, and the 
Higgs mixing angle a are given in this approximation by the relations: 
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as a function of and tan p. The triple Higgs self-couplings of the MSSM can be param- 
eterized in units of M|/f, v = 246 GeV, as. 
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Next-to- leading order (NLO) effects are certainly dominant, though the next-to-next- 
to-leading order (NNLO) ones cannot entirely be neglected (especially in the Higgs mass 
relations). Thus, in the numerical analysis, the complete one-loop and the leading two- loop 
corrections to the MSSM Higgs masses and the triple Higgs self- couplings are included. The 
Higgs masses, widths and self-couplings have been computed using the HDECAY program 
described in Ref. which uses a running 6-mass in evaluating the h ^ bb decay fraction. 
Thus, for consistency, we have evolved the value of rub entering the hbb Yukawa couplings of 
the h —>■ bb decay currents of our processes in the same way. 
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For our analysis, we have considered tan [3 = 3 and 50. For the LHC, high values of tan (3 
produce a signal cross section much larger than the tan f3 = 3 scenario, over almost the entire 
range of m^. However, this enhancement is due to the increase of the down-type quark-Higgs 
coupling, which is proportional to tan (3 itself, and serves only to magnify the dominance of 
the quark box diagrams of Fig. |I]. Unfortunately, these graphs have no dependence on either 
of the two triple Higgs self-couplings entering the gluon-gluon process considered here, i.e., 
Xhhh and Xuhh- Thus, although the cross section is comfortably observable, all sensitivity to 
such vertices is lost. Therefore, the measurement of the triple Higgs self-coupling, Xnhh, is 
only feasible at the LHC for low tan/3 due to the resonant production of the heavy Higgs 
boson (see Fig. 5a of Ref. 0). 

In contrast, the cross section for double Higgs production at the LC is small for large 
tan j3 because there is no heavy Higgs resonance (see Fig. 8 of Ref. [Q]). As soon as it becomes 
kinematically possible to decay the heavy Higgs into a light Higgs pair, the ZZH coupling 
is already too small to generate a sizable cross section. Furthermore, the continuum MSSM 
cross section is suppressed with respect to the SM cross section since the MSSM couplings 
ZZH and ZZh vary with cos(/3 — a) and sin(/5 — a), respectively, with respect to the 
corresponding SM coupling. Notice that in this regime, at a LC, the Xhhh vertex could in 
principle be accessible instead, since Xhhh ^ Xnhh (see Fig. 2 of Ref. 0]) and because of the 
kinematic enhancement induced by ruh <^ mn- Unfortunately, we will see that the size of 
the e'^e~ — ^ hhZ cross section itself is prohibitively small. 

We assume that 6-jets are distinguishable from light-quark and gluon jets and no efficiency 
to tag the four 6-quarks is included in our parton level results. We further neglect considering 
the possibility of the 6-jet charge determination. Also, to simplify the calculations, the Z 
boson appearing in the final state of the LC process is treated as on-shell and no branching 
ratio (BR) is applied to quantify its possible decays. In practice, one may assume that 
it decays leptonically (i.e., Z with ^ = e,fi,T) or hadronically into light-quark 

jets (i.e., Z —>■ qq, with q ^ b), in order to avoid problems with 66-quark combinatorics. 
Furthermore, in the LC analysis, we have not simulated the effects of Initial State Radiation 
(ISR), beamstrahlung or Linac energy spread. Indeed, we expect them to affect signal and 
backgrounds rather similarly, so we can neglect them for the time being. Indeed, since a 
detailed phenomenological study, including both hadronization and detector effects, already 
exists for the case of double Higgs- strahlung in e~^e~ [Q, whose conclusions basically support 
those attained in the theoretical study of Ref. , we limit ourselves here to update the latter 
to the case of the MSSM. 

So far only resonant production gg ^ H hh ^ bbbb has been investigated [1^, with 
full hadronic and detector simulation and considering also the (large) QCD backgrounds, 
and a similar study does not exist for continuum double Higgs production at the LHC. (See 



Ref. [0 for a detailed account of the gg H ~y hh ^ 7766 decay channel.) The event 



simulation has been performed by using a special version of PYTHIA |T8[, in which the 



relevant LO MEs for double Higgs production of Ref. [|TT| have been implemented by M. El 
Kacimi and R. Lafaye. These MEs take into account both continuum and resonant double 
Higgs boson production and their interferences. (The insertion of those for e"'"e~ processes is 
in progress.) The PYTHIA interface to HDECAY has been exploited in order to generate 
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the MSSM Higgs mass spectrum and the relevant Higgs BRs, thus maintaining consistency 
with the parton level approach. As for the LHC detector simulation, the fast simulation 
package was used, with high luminosity (i.e., / Cdt = 100 fb~^) parameters. 

The motivation for our study is twofold. On the one hand, to complement the studies 
of Ref. |T0[ by also considering the continuum production gg —>■ hh bbbb at large tan j3. 



On the other hand, to explore the possibility of further kinematic suppression of the various 
irreducible backgrounds to the resonant channel at small tan j3. 



3 Results 

3.1 The LHC analysis 

In our LHC analysis, following the discussion in Sect. ^, we focus most of our attention on 
the case tan (3 = 3, with = 210 GeV, although other combinations of these two MSSM 
parameters will also be considered. We further set A = — /i = 1 TeV and take all sparticle 
masses (and other SUSY scales) to be as large as 1 TeV. 

3.1.1 gg ^ hh ^ bbbb at parton level 

In our parton level analysis, we identify jets with the partons from which they originate 
(without smearing the momenta) and apply all cuts directly to the partons. We mimic the 
finite coverage of the LHC detectors by imposing a transverse momentum threshold on each 
of the four 6-jets, 

prib) > 30 GeV (4) 
and requiring their pseudorapidity to be 

\r]{b)\ <2.5. (5) 

Also, to allow for their detection as separate objects, we impose an isolation criterium among 
6-jets, 

AR{bb) > 0.4, (6) 

by means of the usual cone variable AR{ij) = Ari{ijy + A(j){ijy, defined in terms of 
relative differences in pseudo-rapidity r]ij and azimuth (pij of the i-th and j-th 6-jets. 

As preliminary and very basic selection cuts (also to help the stability of the numerical 
integration), we have required that the invariant mass of the entire 46-system is at least 
twice the mass of the lightest MSSM Higgs boson (apart from mass resolution and gluon 
emission effects), e.g., 

m{bbbb) > 2m/, - 40 GeV, (7) 
and that exactly two /i-resonances are reconstructed, such that 

|m(66) -m,,| < 20 GeV. (8) 
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Figure 2: Distributions in transverse momentum of the four pj'-ordered 6-jets in gg ^ hh 
bbbb and in the QCD background, after the cuts at the LHC, for tan/3 = 3, rrih = 104 

GeV and rriH — 220 GeV. Normahzation is to unity. 

In doing so, we imphcitly assume that the h scalar boson has aheady been discovered and 
its mass measured through some other channel, as we have already intimated in the Intro- 
duction. 

After the above cuts have been implemented, we have found that the two 46-backgrounds 
proceeding through EW interactions are negligible compared to the pure QCD process. In 
fact, the constraints described in eqs. produce the strongest suppression, almost 

completely washing out the relatively enhancing effects that the cuts in (|^)-(|^) have on 
the EW components of the backgrounds with respect to the pure QCD one, owning to the 
intermediate production of massive Z bosons in the former. In the end, the production 
rates of the three subprocesses scale approximately as their coupling strengths: i.e., C(ag) : 
O^a^alj^) : 0{a'^^). Therefore, in the reminder of our analysis, we will neglect EW effects, 
as they represent not more than a 10% correction to the QCD rates, which are in turn 
affected by much larger QCD A'-factors. As for the pure QCD background itself, it hugely 
overwhelms the double Higgs signal at this stage. The cross section of the former is about 
7.85 pb, whereas that of the latter is approximately 0.16 pb. 

To appreciate the dominance of the rrih cuts, one may refer to Fig. ^ where the distri- 
butions in transverse momentum of the four pj^-ordered 6-quarks (such that pribi) > ... > 
Pribi)) of both signal and QCD background are shown. Having asked the four 6-jets of the 
background to closely emulate the gg hh bbbb kinematics, it is not surprising to see 
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a 'degeneracy' in the shape of all spectra. Clearly, no further background suppression can 
be gained by increasing the prib) cuts. The same can be said for ri{b) and AR{bb). Others 
quantities ought to be exploited. 




Figure 3: Distributions in minimum relative angle (in radians) in the 46-system rest frame 
between two 6-jets reconstructing rrih in gg ^ hh —>■ bbbb and in the QCD background, 
after the cuts (|)-(|) at the LHC, for tan/3 = 3, = 104 GeV and tuh ^ 220 GeV. 
Normalization is to unity. 

In Fig. 0, we present the signal and QCD background distributions in the minimum angle 
formed between the two 6-quarks coming from the 'same Higgs' (i.e., those fulfilling the cuts 
in (^) in the 46-system rest frame (the plot is rater similar for the maximum angle, thus 
also on average). There, one can see a strong tendency of the two 26-pairs produced in the 
Higgs decays to lie back-to-back, reflecting the 2-^2 intermediate dynamics of Higgs pair 
production via gg — » hh. Missing such kinematically constrained virtual state, the QCD 
background shows a much larger angular spread towards small ^min(&&) values, eventually 
tamed by the isolation cut (H). 

The somewhat peculiar shape of the signal distribution is due to destructive interference. 
Recall that the signal contains not only diagrams proceeding via a heavy Higgs resonance (the 
upper-left hand graph of Fig. |I]), which results in the large peak in Fig. |^, but also contains 
a continuum contribution mediated by box graphs (the upper-right hand graph of Fig. |l]). 
These two contributions destructively interfere leading to the depletion of events between 
the large back-to-back peak and the small remaining 'bump' of the continuum contribution 
as seen in Fig. |^. 

In the end, a good criterium to enhance the signal-to-background ratio {S/ B) is to require, 
e.g., 6{bb) > 2.4 radians, i.e., a separation between the 26-jets reconstructing the lightest 
Higgs boson mass of about 140 degrees in angle. (Incidentally, we also have investigated the 
angle that each of these 26-pairs form with the beam axis, but found no significant difference 
between signal and QCD background). 
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Figure 4: Distributions in thrust in the rest frame of the 46-system in gg ^ hh bbbb and 
in the QCD background, after the cuts ®-(|D at the LHC, for tan /? = 3, rrih = 104 GeV 
and TJiH — 220 GeV. Normahzation is to unity. 



An additional consequence that one should expect from the presence of two intermediate 
massive objects in gg ^ hh ^ bbbb events is the spherical appearance of the jets in the final 
state, in contrast to the usual planar behavior of the infrared QCD interactions. These phe- 
nomena can be appreciated in Fig. ^. Notice there the strong tendency of the background to 
yield high thrust configurations, again controlled by the separation cuts when T approaches 
unity. On the contrary, the average value of the thrust in the signal is much lower, being the 
effect of accidental pairings of 'wrong' 26-pairs (the shoulder at high thrust values) marginal. 
An effective selection cut seems to be, e.g., T < 0.85. 

Furthermore, if the heavy Higgs mass is sufficiently well measured at the LHC then one 
can exploit the large fraction of 46-events which peak at rriH in the signal, as dictated 
by the H ^ hh decay, improving the signal-to-background ratio. This peak at ttlh can be 
clearly seen in the left hand plot of Fig. ^, where it dominates the QCD background, even for 
bins 13 GeV wide. In fact, not only could the QCD background be considerably suppressed 
but also those contributions to gg hh not proceeding through an intermediate H state 
should be removed, this greatly enhancing the sensitivity of the signal process to the Xnhh 
coupling. This can be seen in the right hand plot of Fig. ^ where the signal is shown on a 
logarithmic scale. The continuum contribution due to the box graphs (and its destructive 
interference with the heavy Higgs decay contribution) is now evident although one should 
note that it is considerably suppressed compared to the peak at mn- 

Now, if a less than 10% mass resolution can be achieved on the light and heavy Higgs 
masses, then one can tighten cut (§) to \m{bb) — mh\ < 10 GeV and introduce the additional 
cut \m{bbbb) — mff | < 20 GeV. These cuts taken together with those in 6{bb) and T already 
suggested, reduce the QCD background to the same level as the signal. In fact, we have 
found that the cross section of the background drops to approximately 174 fb whereas that 
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Figure 5: Distributions in invariant mass of the 46-system in gg ^ hh ^ bbbb and in the 
QCD background, after the cuts (Q)-® at the LHC, for tan/? = 3, = 104 GeV and 
rriH — 220 GeV. Normahzation is to unity. The left hand plot shows both the signal (solid 
curve) and the QCD background (dashed curve), distributed in 5 GeV bins. The same signal 
is also shown as a histogram for a more experimentally realistic binning of 13 GeV. The right 
hand plot also shows the signal (collected in 5 GeV wide bins) on a logarithmic scale. Here 
the structure of the continuum contribution (and its destructive interference with the heavy 
Higgs decay contribution) can be seen. 

of the signal remains as large as 126 fb, this yielding a very high statistical significance at 
high luminosity. Even for less optimistic mass resolutions the signal-to-background ratio 
is still significantly large. For example, selecting events with \m{bb) — mh\ < 20 GeV and 
\m{bbbb) — mnl < 40 GeV, the corresponding numbers are approximately 102 fb for the 
signal and 453 fb for the background. Notice that the signal actually decreases as these 
Higgs mass windows are made larger. This is due to our insistence that exactly two bb pairs 
should reconstruct the light Higgs mass. As the light Higgs mass window is enlarged from 
rrih ± 10 GeV to rrih ± 20 GeV, it becomes more likely that accidental pairings reconstruct 
the light Higgs boson. Since one is then unable to unambiguously assign the b quarks to the 
light Higgs bosons, the event is rejected and the signal drops. 

Although we have discussed here an ideal situation which is difficult to match with more 
sophisticated hadronic and detector simulations, it still demonstrates that the measurement 
of the Xuhh coupling could be well within the potential of the LHC, at least for our particular 
choice of MSSM parameters. Comforted by such a conclusion, we now move on to more 
realistic studies. 

3.1.2 gg ^ hh ^ bbbb at the LHC experiments 

Although the LHC experiments will be the first where one can attempt to measure the Higgs 
self-couplings, the analysis is very challenging because of the smallness of the production cross 
sections. Even in the most favorable cases, the production rate is never larger than a few 
picobarns, already including one-loop QCD corrections, as computed in Ref. [|ri|] . The cross 
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sections at this accuracy are given in Tab. ^ for the resonant process (case 1 with mn = 220 
GeV) as well as three non resonant scenarios: one at the same tan (3 but with the H hh 
decay channel closed (case 2), a second at very large tan/5 and no visible resonance (case 3) 
and, finally, the SM option (case 4, where rrih identifies with the mass of the standard Higgs 
state). 



case 


model 


tan P 


mh (GeV) 


A (TeV) 


/i (TeV) 


a(fb) 


dominant mode 


1 


MSSM 


3 


104 


+ 1 


-1 


2000 


gg H ^ hh 


2 


MSSM 


3 


100 


+1 


-1 


20 


gg hh 


3 


MSSM 


50 


105 


+ 1 


+ 1 


5000 


gg hh 


4 


SM 




105 






40 


gg ^ hh 



Table 1: Cross sections for double Higgs production hh at the LHC via gluon-gluon fusion 
at NLO accuracy, for three possible configurations of the MSSM and in the SM as well. 

3.1.3 LHC trigger acceptance 

For 46-final states, possible LHC triggers are high pt electron/muons and jets. As an exam- 
ple, the foreseen ATLAS level 1 trigger thresholds on px and acceptance for a 46-selection 
(with the four 6-jets reconstructed in the detector) are given in Tab. ^, assuming the LHC 
to be running at high luminosity. 



trigger type: 


1 e 


1 fi 


2 /i 


1 jet 


3 jets 


4 jets 


total 


Pt in GeV 


30 


20 


10 


290 


130 


90 




case 1, e{bbbb) in % 


0.01 


0.01 


0.4 


0.08 


0.08 


0.05 


0.53 


case 2 


< 0.01 


< 0.01 


2.1 


2.9 


3.8 


4.2 


8.8 


case 3 


< 0.01 


< 0.01 


2.2 


2.7 


3.8 


4.1 


8.7 


case 4 


< 0.01 


< 0.01 


2.0 


2.5 


3.3 


3.6 


7.8 



Table 2: Kinematical acceptance of the ATLAS detector to trigger four 6-jets (including 
detector acceptance) at high luminosity. 

The overall trigger acceptance is at best 8-9%, for cases 2,3,4. The very low efficiency for 
case 1 is clearly a consequence of the small value of the difference mu — 2mh, translating into 
a softer pT{b) spectrum with respect to the other cases (compare the left-hand with the right- 
hand side of Fig. One can further see in the left-hand plot of Fig. ^ that the bulk of the 
signal lies below the lowest prib) threshold of Tab. || (i.e., 90 GeV), so that adopting smaller 
trigger thresholds could result in a dramatic enhancement of our efficiency. Of course, this 
would also substantially increase the low transverse momentum QCD background, as we can 
see in the parton level analysis of Fig. |^. 

For example, by lowering the thresholds to 180, 80 and 50 GeV for 1, 3 and 4 jets, 
respectively (compare to Tab. H), the overall trigger acceptance on the signal goes up to 



77 





50 100 ISO 200 250 300 50 100 150 200 250 300 

b-jetslransverseenergy E^CGeV) b-jetstransverseenergy E^(GeV) 

Figure 6: Reconstructed transverse energy/momentum for 6-jets in gg ^ hh ^ bbbb events 
of case 1 (left plot) and 6-jets in gg ^ hh bbbb events of case 2 (right plot) with ATLAS 
fast simulation at high luminosity. Normalization is arbitrary. 

1.8%, i.e., by almost a factor of 4. Meanwhile though, the ATLAS level-1 jet trigger rates 
increase by a factor of 10 ||19|. Anyhow, even for our poor default value of e{bbbb) in Tab. ^ 
we will see that case 1 still yields a reasonable number of events in the end. Optimizations 
of the 6-jet transverse momentum thresholds are in progress [Q. 

3.1.4 LHC events selection for gg ^ hh ^ bbbb 

Jets are reconstructed merging tracks inside AR{bb) = 0.4. Only jets with transverse en- 
ergy/momentum greater than 30 GeV and with \ri{b)\ < 2.5 are kept. (Thus, the same 
cuts as in the parton level analysis, now applied instead to jets.) The effect from pile up is 
included in the resolution. A jet energy correction is then applied. 

The invariant masses of each jet pair can then be computed. Assuming that the lightest 
Higgs boson mass is known, events with m{bb) sufficiently close to rrih can efficiently be 
selected, see Fig. |^. Another cut on the AR{bb, bb) between pairs of 6-jets can also be 
applied to reduce the intrinsic combinatorial background, since the latter concentrates at 
large AR{bb, bb) values, see Fig. |^. 

For case 1, as already discussed, we can further impose that the invariant mass of the 
four 6-jets should be the heavy Higgs mass, tuh, in order to select the H ^ hh resonance, 
as confirmed by Fig. ^. In the other three cases, where the H ^ hh splitting is no longer 
dominant (MSSM) or non-existent (SM), one can still insist that the 46-jet invariant mass 
should be higher than two times the lightest Higgs mass, see Fig. |1^ and recall eq. (|^. 
Finally, following Fig. |TT|, by constraining the 6-jets pairs four-momenta around the known 
light Higgs mass value, m/j, one can further reject the intrinsic background by means of the 
m{bbbb) spectrum. 
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Figure 7: Reconstructed invariant mass distribution of 26-jet pairs in continuum gg 
hh — > bbbb events (case 2) with the fast simulation at high luminosity. Normalization is 
arbitrary. (Results of a Gaussian fit are also given.) 
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Figure 8: Reconstructed AR{bb, bb) between 26-jet systems from h ^ bb decays in contimmm 
gg ^ hh ^ bbbb events (case 2) with the fast simulation at high luminosity. The dashed 
histogram shows the same distribution for all pairs of jets. Normalization is arbitrary. 
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Figure 9: Reconstructed 46-jet invariant mass for 6-jets coming from the hh pair in gg 
hh bbbb events (case 1) with the fast simulation at high luminosity. The dashed histogram 
shows the same distribution for all groups of four jets. Normalization is arbitrary. (Results 
of a Gaussian fit to the first spectrum are also given.) 

3.1.5 LHC 6-tagging in gg ^ hh ^ bbbb 

The 6-tagging efficiency at high luminosity is set to 50%, with pj- dependent correction 
factors for jets rejection. An average rejection of 10 for c-jets and 100 for light-jets can be 
expected. We then studied the effect on the selection efficiency of requiring from one to four 
6-tags, although it is clear that, according to the parton level studies, the huge background 
rate demands four 6-tags, leading to a 6% tagging efficiency overall. 

3.1.6 Event rates at the LHC 

Taking into account all the efficiencies described above, and using the NLO normalization of 
Tab. m, one can extract the number of expected events per year at the LHC at high luminosity 
given in Tab. The selection cuts enforced here are the following. For a start, we have kept 
configurations where \m{bb) — mfi\ < 30 GeV (cases 1,3,4) or \m{bb) — mfi\ < 20 GeV (case 2) 
and AR{bb, bb) < 2.5 (all four cases). (If more than two m/^'s are reconstructed, the best two 
26-pairs are selected according to the minimum value of = [mh—m{bb)]'^ + [mh—m'{bb)'^].) 
Then, a cut on m{bbbb) is applied: in presence of the H ^ hh resonance (case 1) we have 
kept events within an mn mass window of ±2cr (about 82% of the total number survive); 
otherwise (cases 2,3,4) we have adjusted the m{bbbb) > 2mh cut so to keep 90% of the sample. 
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Figure 10: Reconstructed 46-jet invariant mass for 6-jets coming from the hh pair in gg 
hh bbbb events (case 4) with the fast simulation at high luminosity. The dashed histogram 
shows the same distribution for all groups of four jets. Normalization is arbitrary. 

In the end, one finds the numbers in Tab. ^ that are encouraging indeed. 





case 1 


2 


3 


4 


cr in fb 


2000 


20 


5000 


40 


trigger threshold acceptance 


0.53% 


8.8% 


8.7% 


7.8% 


mass windows 


60% 


50% 


40% 


40% 


46-tagging 


6% 


6% 


6% 


6% 


events/year (no tagging) 


636 


88 


17400 


125 


events/year (four 6-tags) 


38 


5.3 


1044 


7.5 



Table 3: Total rates for gg —>■ hh ^ bbbb, after all efficiencies have been included and 
selection cuts (§)-(^ enforced at hadron level, with 100 fb~^ per year of luminosity. 

In conclusion then, looking at the results in Tab. |^ and bearing in mind the potential seen 
in reducing the pure QCD background via gg 0{a^) bbbb (see Figs. |]-|), one should be 
confident in the LHC having the potential to measure the Xnhh coupling in resonant H ^ hh 
events (case 1). To give more substance to such a claim, we have now initiated background 
studies at hadron and detector level, following the guidelines obtained by the parton level 
analysis 0. As for other configurations of the MSSM (such as case 2) or in the SM (case 4), 
the expectations are more pessimistic. Case 3 deserves further attention. In fact, notice the 
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Figure 11: Reconstructed 46-jet invariant mass for 6-jets coming from the hh pair in gg — > 
hh — >■ bbbb events (case 4) with the fast simulation at high luminosity. Here, the energy of 
the jet pairs is recalculated using the rrih constraint. The dashed histogram shows the same 
distribution for all groups of four jets. Normalization is arbitrary. 

large number of events surviving and recall what mentioned in the Introduction concerning 
the potential of the non-resonant gg ^ hh bbbb process as a discovery channel of the light 
Higgs boson of the MSSM in the large tan f3 region at moderate rriA values, a corner of the 
parameter space where the h coverage is given only by SM-like production/decay modes, 
thus not allowing one to access information on the MSSM parameters. Results on this topic 
too will be presented in Ref. 

3.2 The LC analysis 

Here, we closely follow the selection procedure advocated in Ref. . In order to resolve the 
four 6-jets as four separate systems inside the LC detector region, we impose the following 
cuts. First, that the energy of each 6-jet is above a minimum threshold, 

E{b) > 10 GeV. (9) 

Second, that any 6-jet is isolated from all others, by requiring a minimum angular separation, 

cos ^(6, 6) < 0.95. (10) 
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Similarly to the hadronic analysis, one can optimize S/B by imposing the constraints [Q, 



m{hhhh) > 2mh - 10 GeV, (11) 

|m(66) -m,,| < 5 GeV, (12) 

on exactly two combinations of 26-jets. Here, note that the mass resolution adopted for the 
quark systems is significantly better than in the LHC case, due to the cleanliness of the e"'"e~ 
environment and the expected performance of the LC detectors in jet momentum and angle 
reconstruction E^. Thus, given such high mass resolution power from the LC detection 



apparatus, one may further discriminate between h and Z mass peaks by requiring that 
none of the 26-jet pairs falls around m^. 



\m{bb) — mz\ > 5 GeV. 



(13) 



Moreover, in the double Higgs-strahlung process e~^e~ — hhZ, the four 6-quarks are pro- 
duced centrally, whereas this is generally not the case for the background (see the discussion 
in Ref. f^). This can be exploited by enforcing 



\ cos 9{bb,bbb,bbbb)\ < 0.75, 
where 6{bb, bbb, bbbb) are the polar angles of all two-, three- and four-jet systems. 



(14) 
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Figure 12: Cross sections in femtobarns for the e+e^ hhZ signal in the h bbbb decay 
channel, at a LC with 500 GeV as CM energy, as a function of rrih for tan/? = 3 and 
50. Our acceptance cuts in energy and separation of the four 6-quarks (PD-([TO|) have been 
implemented. No beam polarization is included. 



Fig. |I2| shows the production and decay rates of the signal process, e e — > hhZ bbbbZ, 
as obtained at the partonic level, after the cuts (|9|)- (p!o|) have been implemented. The MSSM 
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setup here includes some mixing, having adopted A = 2.4 TeV and /i = 1 TeV, at both 
tan/3 = 3 and 50. Notice the onset of the H ^ hh ^ bbbb decay sequence in the Higgs- 
strahlung process e~^e~ — > HZ at low tan/5. The same does not occur for large values, as 
previously explained. The impact of the above jet selection cuts on the signal is marginal, 
as the 6-quarks are here naturally isolated and energetic, being the decay products of heavy 
objects. In fact, the rates in Fig. |12| would only be 10-20% higher if all the 46-quark phase 
space was allowed (the suppression being larger for smaller Higgs masses). At the height 
of the resonant peak around rrih ~ 104 GeV at tan /3 = 3, the signal rate is not large but 
observable, yielding more than one event every 1 fb~^ of data. For a high luminosity 500 
GeV TESLA design p2|], this would correspond to more than 300 events per year. Given 



the very high efficiency expected in tagging 6-quark jets, estimated at 90% for each pairs of 
heavy quarks one should expect a strong sensitivity to the triple Higgs self-coupling. 
The situation at large tan (3 is much more difficult instead, being the production rates smaller 
by about a factor of 10. 



In the left-hand side of Fig. we present the EW background, after the constraints in 
(P|)-(PI]|) have been enforced, in the form of the four dominant EW sub-processes. These four 
channels are the following. 

1. e~^e^ —>■ ZZZ —>■ bbbbZ, first from the left in the second row of topologies in Fig. 3 of 
Ref. That is, triple Z production with no Higgs boson involved. 

2. e+e^ h/HZZ —>■ bbbbZ, first(first) from the left(right) in the fifth(fourth) row of 
topologies in Fig. 2 of Ref. ||^ (also including the diagrams in which the on-shell Z is 
connected to the electron-positron line). That is, single Higgs-strahlung production in 
association with an additional Z, with the Higgs decaying to bb. The cross sections of 
these two channels are obviously identical. 

3. e~^e^ — > h/HZ Z*Z*Z bbbbZ, first from the right in the third row of topologies 
in Fig. 2 of Ref. . That is, single Higgs-strahlung production with the Higgs decaying 
to bbbb via two off-shell Z* bosons. 

4. e~^e~ Zh/H bbZ*Z bbbbZ, first(first) from the right(left) in the first(second) 
row of topologies in Fig. 2 of Ref. ^ . That is, two single Higgs-strahlung production 
channels with the Higgs decaying to bbZ via one off-shell Z* boson. Also the cross 
sections of these two channels are identical to each other, as in 2. 

The Oi^a^al^) EW/QCD background is dominated by e~^e~ ZZ production with one 
of the two Z bosons decaying hadronically into four 6-jets. This subprocess corresponds to 
the topology in the middle of the first row of diagrams in Fig. 4 of Ref. P]. Notice that 
Higgs graphs are involved in this process as well (bottom-right topology in the mentioned 
figure of [^). These correspond to single Higgs-strahlung production with the Higgs scalar 
subsequently decaying into bbbb via an off-shell gluon. Their contribution is not entirely 
negligible, owing to the large ZH production rates, as can be seen in the right-hand side of 
Fig. O. The interferences among non-Higgs and Higgs terms are always negligible. 
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Figure 13: Cross sections in femtobarns for the dominant components of the EW (left) and 
EW/QCD (right) background to the e~^e~ hhZ signal in the h bbbb decay channel, 
at a LC with 500 GeV as CM energy, as a function of for tan /3 = 3 (top) and 50 
(bottom). Our acceptance cuts in energy and separation of the four 6-quarks (P)-(p!0|) have 
been implemented. No beam polarization is included. 

In performing the signal-to-background analysis, we have chosen two representative points 
only, identified by the two following combinations: (i) tan (3 = 3 and = 210 GeV 
(yielding rrih ~ 104 GeV and mn ~ 220 GeV); (ii) tan/3 = 50 and = 130 GeV (yielding 
rrih ~ 120 GeV and mn ~ 130 GeV). These correspond to the two asterisks in Fig. ^ 
that is, the maxima of the signal cross sections at both tan f3 values. The first corresponds 
to resonant H hh production, whereas the latter to the continuum case. If we enforce 
the constraints of eq. (|TT])-(|np, the suppression of both EW and EW/QCD is enormous, so 
that the corresponding cross sections are of C(10~^) fb, while the signal rates only decrease 
by a factor of four at most. This is the same situation that was seen for the SM case in 
Ref. 0. Indeed, in the end it is just a matter of how many signal events survive, the sum 
of the backgrounds representing no more than a 10% correction (see Fig. 11 of Ref. |^). 
For example, after 500 fb~^ of data collected, one is left with 156 and 15 events for case (i) 
and (ii), respectively. However, these numbers do not yet include 6-tagging efficiency and Z 
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decay rates. 



4 Summary 

To summarize, the 'double Higgs production' subgroup has contributed to the activity of 
the Higgs WG by assessing the feasibihty of measurements of triple Higgs self-couplings at 
future TeV colliders. The machines considered were the LHC at CERN (14 TeV) and a 
future LC running at 500 GeV. In both high luminosity setup was assumed, given 

the smallness of the double Higgs production cross sections. In particular, the H hh 
resonant enhancement was the main focus of our studies, involving the lightest, h, and the 
heaviest, H, of the neutral Higgs bosons of the MSSM, in the kinematic regime ttlh > 2mh- 
This dynamics can for example occur in the following reactions: gg — > hh in the hadronic 
case and e'^e~ — >■ hhZ in the leptonic one, but only at low tan/?. These two processes 
proceed via intermediate stages of the form gg H and e~^e~ HZ, respectively, followed 
by the decay H hh. Thus, they in principle allow one to determine the strength of the 
Hhh vertex involved, Xuhh-, in turn constraining the form of the MSSM Higgs potential itself. 
The signature considered was hh bbbb, as the h ^ bb decay rate is always dominant. 

We have found that several kinematic cuts can be exploited in order to enhance the 
signal-to-background rate to level of high significance, particularly at the e~^e~ machine. At 
the pp accelerator, in fact, the selection of the signal is made much harder by the presence 
of an enormous background in 46 final states due to pure QCD. In parton level studies, 
based on the exact calculation of LO scattering amplitudes of both signals and backgrounds 
(without any showering and hadronization effects but with detector acceptances), we have 
found very encouraging results. At a LC, the double Higgs signal can be studied in an 
essentially background free environment. At the LHC, the signal and the QCD background 
are in the end at the same level with detectable but not very large cross sections. 

Earlier full simulations performed for the e~^e~ case had already indicated that a more so- 
phisticated treatment of both signal and backgrounds, including fragmentation/hadronization 
and full detector effects, should not spoil the results seen at the parton level. For the LHC, 
our preliminary studies of gg ^ H ^ hh ^ bbbb in presence of the gg hh bbbb 
continuum (and relative interferences) also point to the feasibility of the signal selection, 
after realistic detector simulation and event reconstruction. As for double h production in 
the continuum, although not very useful for Higgs self-coupling measurements, this seems 
a promising channel, if not to discover the lightest MSSM Higgs boson certainly to study 
its properties and those of the Higgs sector in general (because of the large production and 
decay rates at high tan f] and its sensitivity to such a parameter), as shown from novel simu- 
lations also presented in this study. (The discovery potential of this mode will eventually be 
addressed in Ref. 0.) Despite lacking a full background analysis in the LHC case, we have 
no reason to believe that a comparable degree of suppression of background events seen at 
parton level cannot be achieved also at hadron level. Progress in this respect is currently 
being made 0. 
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Programs and Tools for Higgs Bosons 

E. Boos, A. Djouadi, N. Ghodbane, S. Heinemeyer, 
V. Ilyin, J. Kalinowski, J.L. Kneur and M. Spira 

Abstract 

The search strategies for Higgs bosons at LEP, Tevatron, LHC and future e'^e~ 
hnear coUiders (LC) and muon coUiders exploit various Higgs boson production and 
decay channels. The strategies depend not only on the experimental setup [e.g. hadron 
versus lepton colliders] but also on the theoretical scenarii, for instance the Standard 
Model (SM) or some of its extensions such as the Minimal Supersymmetric Standard 
Model (MSSM). It is of vital importance to have the most reliable predictions for the 
Higgs properties, branching ratios and production cross sections. 

There exist several programs and packages which determine the properties of Higgs 
particles, their decays modes and production mechanisms at various colliders. These 
programs are in general independent, have different inputs and treat different aspects 
of the Higgs profile. During this workshop, many discussions have been made and 
some work has been done on how to update these various programs to include the 
latest theoretical developments, and how to link some of them. 

This report summarizes the work which has been performed in this context. 

1 HDECAY 

The program HDECAY can be used to calculate Higgs boson partial decay widths and 
branching ratios within the SM and the MSSM and includes: 

• All decay channels that are kinematically allowed and which have branching ratios 
larger than 10~^, y compris the loop mediated, the three body decay modes and in the 
MSSM the cascade and the supersymmetric decay channels 0. 

• In the MSSM, the complete radiative corrections in the effective potential approach 
with full mixing in the stop/sbottom sectors; it uses the renormalization group im- 
proved values of the Higgs masses and couplings and the relevant next-to-leading- 
order corrections are implemented 0. 

• All relevant higher-order QCD corrections to the decays into quark pairs and to the 
loop mediated decays into gluons and photons are incorporated in a complete form ^ ; 
the small leading electroweak corrections are also included. 

• Double off-shell decays of the CP-even Higgs bosons [SM Higgs and the h, H bosons 
of the MSSM] into massive gauge bosons which then decay into four massless fermions, 
and all important below-threshold three-body decays [decays into one real and virtual 
gauge bosons, cascade decays into a Higgs and a virtual gauge boson, decays into a 
real and virtual heavy top quark, etc,..] 
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• In the MSSM, all the decays into SUSY particles [neutralinos, charginos, sleptons 
and squarks including mixing in the stop, sbottom and stau sectors] when they are 
kinematically allowed 0. 

• In the MSSM, the SUSY particles are also included in the loop mediated 77 and gg 
decay channels, with the leading parts of the QCD corrections incorporated |^. 

The source code of the program, hdecay.f written in FORTRAN, has been tested on 
computers running under different operating systems. It is self-contained and all the nec- 
essary subroutines [e.g. for integration] are included. The program provides a very flexible 
and convenient usage, fitting to all options of phenomenological relevance. The program 
is lengthy [more than 6000 lines] but rather fast, especially if some options [as decays into 
double off-shell gauge bosons] are switched off. 

The basic input parameters, fermion and gauge boson masses and their total widths, 
coupling constants and, in the MSSM, soft SUSY-breaking parameters can be chosen from 
an input file hdecay.in. In this file several flags allow switching on/off or changing some 
options [e.g. choosing a particular Higgs boson, including/excluding the multi-body or 
SUSY decays, or including/excluding specific higher-order QCD corrections]. 

The results for the many decay branching ratios and the total decay widths are written 
into output files br.Xi [with X = H^, h, H, A and i = 1,...] with headers indicating the 
various processes and giving some of the parameters. 

Since the release of the original version of the program several bugs have been fixed and 
a number of improvements and new theoretical calculations have been implemented. During 
this workshop, the following points have been included: 

• Link to the FeynHiggsFast routine which gives the masses and couplings of the MSSM 
up to two-loop order in the diagrammatic approach |^. 

• Link to the SUSPECT routine for the Renormalisation Group evolution and for the 
proper electroweak symmetry breaking in the minimal Supergravity model [^. 

• Implementation of Higgs boson decays to a gravitino and neutralino or chargino in 
gauge-mediated SUSY breaking models |]10[ . 

• Inclusion of gluino loops in Higgs boson decays to qq pairs [pT| . 

• Determination and inclusion of the RG improved two-loop contributions to the MSSM 
Higgs boson self-interactions. 

In addition, the inclusion of the [possibly large] QCD corrections for the MSSM Higgs 



boson decays into squark pairs |12| has started. 

The log-book of all modifications and the most recent version of the program can be 
found on the web page http://www.desy.de/~spira/prog. 
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2 Programs for Higgs production 



Several programs for Higgs boson production at hadron colliders in the context of the SM 
and the MSSM, including the next-to-leading order (NLO) QCD corrections, are available 
at the web page: http : //www . desy . de/~spira. The purpose of these programs, and some 
improvements made during this Workshop, are summarized below. For the physics context, 
see the contribution in Section 5 of these proceedings. 

HIGLU calculates the total cross sections for Higgs production in the gluon-fusion mech- 
anism, gg — > Higgs, including the NLO QCD corrections in the SM, MSSM and in a general 
two-Higgs doublet model [by initializing the Yukawa couplings to quarks]. It includes both 
top and bottom quark loops which generate the Higgs couplings to gluons. Moreover the 
program calculates the decay width of Higgs bosons into gluons at NLO. 

V2HV calculates the LO and NLO cross sections for the production in the Higgs-strahlung 
mechanism, qq ^ V + ^ where V — W/Z and $ is a CP-even Higgs boson. The QCD 
corrections are those of the Drell-Yan process; see Section 5. 

VV2H calculates the LO and NLO cross sections for the production in the weak vector 
boson fusion mechanism, qq V*V* ^ qq^ where $ is a CP-cvcn Higgs boson. The QCD 
corrections are included in the structure function approach; see Section 5. 

HQQ calculates the LO cross sections for the production of neutral Higgs bosons in asso- 
ciation with heavy quarks, gg/qq QQ+ Higgs. The NLO QCD corrections are not yet 
completely available and are not included. 

HPAIR calculates the LO and NLO cross sections for the production of pairs of neutral 
Higgs bosons in the the gluou gluon fusion mechanism, gg — > $i$25 or in the Drell-Yan like 
process, qq — > $i$2- The NLO corrections are included only in the heavy top quark limit 
for the gg process. 

The source programs are written in FORTRAN and have been tested on computers run- 
ning under different operating systems. In most cases, the various relevant input parameters 
can be chosen from an input file including a flag specifying the model. 

Since the first release of these programs, the following improvements have been made 
[some of them during this Workshop] : 

• A link to different subroutines calculating the MSSM Higgs boson masses and couplings 
has been installed for all the programs. 

• The contribution of squark loops has been included in HIGLU. 

• The SUSY-QCD corrections have been included in V2HV and VV2H. 

• The contribution of initial 6-quark densities has been included in HQQ. 

• The new version of HDECAY for the neutral Higgs boson total decay widths has been 
included in HPAIR. 
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3 FeynHiggsFast 



In this section|3 we present the Fortran code FeynHiggsFast. Starting from low energy 
MSSM parameters [rrit the top quark mass, tan/3 the ratio of the vev's of the two Higgs 
doublets, the pseudoscalar Higgs mass Ma, the soft SUSY breaking scalar masses M^^, M^^, 
the trilinear coupling At and the higgsino mass parameter n] , FeynHiggsFast calculates the 
masses of the neutral CP-even Higgs bosons, Mh and Mh, as well as the corresponding 
mixing angle a, at the two-loop level In addition the mass of the charged Higgs boson, 
Mh±, is evaluated at the one-loop level. The p-parameter, which allows for constraints in 
the scalar fermion sector of the MSSM, is evaluated up to 0{aas), taking into account the 
gluon exchange contribution at the two-loop level [|T^. 

FeynHiggsFast is based on a compact analytical approximation formula, containing at 
the two-loop level the leading corrections of 0{aas) obtained in the Feynman-diagrammatic 
approach and of (9((j|^m^) obtained with renormalization group (RG) methods Con- 
trary to the full result in the FD approach [§] which has been incorporated into the FOR- 



TRAN code FeynHiggs the approximation formula is much shorter. Thus, the program 
FeynHiggsFast is about 3 x 10^ times faster than FeynHiggs, while the agreement between 
the two codes is better than 2 GeV for the CP-even Higgs bosons masses in most parts of 
the MSSM parameter space. 

The complete program FeynHiggsFast consists of about 1300 lines FORTRAN code. 
The executable file fills about 65 KB disk space. The calculation for one set of parameters, 
including the Ap constraint, takes about 2 x 10~^ seconds on a Sigma station [Alpha pro- 
cessor, 600 MHz processing speed, 512 MB RAM]. The program can be obtained from the 
FeynHiggs homepage: http://www-itp.physik.uni-karlsruhe.de/feynhiggs where 
the code itself is available, together with a short instruction, information about bug fixes, 
etc... 

FeynHiggsFast consists of a front-end, program FeynHiggsFast, and the main part 
where the calculation is performed, starting with subroutine f eynhiggsf astsub. The 
front-end can be manipulated by the user at will, whereas the main part should not be 
changed. In this way FeynHiggsFast can be accommodated as a subroutine to existing 
programs, thus providing an extreme fast evaluation for the masses and mixing angles in the 
MSSM Higgs sector. As discussed previously, this has already been successfully performed 
for the program HDECAY during this workshop. 

FeynHiggsFast asks for the low energy SUSY parameter, listed in Table 1. Concerning 
the stop sector, the user has the option to enter either the physical parameters, i.e. the 
masses and the mixing angle (m^^ , m^^ and sin 6^) or the unphysical soft SUSY breaking 
scalar mass parameters M^^, M^^ and the mixing parameter M/"^ = mt{At — fi cot (3). From 
these input parameters FeynHiggsFast calculates the masses and the mixing angle of the 
MSSM neutral CP-even Higgs bosons, as well as the mass of the charged Higgs boson and 
the p parameter. 
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Table 1: The meaning of the different MSSM variables to be entered into FeynHiggsFast. 

4 SUSPECT 

The fortran codeQ SUSPECT calculates essentially the masses and some of the couplings 
of the SUSY and Higgs particles within the framework of the MSSM. It includes several 
specific options whose purpose is, hopefully, to gain more flexibility with the generally non- 
trivial Lagrangian-to-physical parameter relationship in the MSSM. In particular, besides 
the now widespread procedure of evolving the soft parameters from some universal "min- 
imal SUGRA" high energy initial values down to obtain a corresponding low-energy spec- 
trum, SUSPECT can also treat almost arbitrary non-universal departures from this SUGRA 
model. The latest version 1.2 is a subroutine, so that it can be easily interfaced with any 
other FORTRAN codes, as will be described below. It also includes some new useful tools 
like, for instance, the possibility of evolving the parameters " bottom-up", the possibility of 
choosing as input some of the parameters that are usually obtained as output, etc. 

The latest version of the program consists of three parts: the subroutine suspectl2.f , 
suspectl2-call . f an example of calling routine and suspectl2.in a typical example of 
input file. To interface SUSPECTl . 2 properly with your own main code, the easiest way is 
first to run the example code suspectl2-call . f . Once familiar with the calling procedure, 
you may simply implement in your calling code a few appropriate command lines stripped 
from the example file, that you can adapt to your purpose. 

The core of the SUSPECT algorithm is conveniently separated into three different tasks, 
that are indeed conceptually -and technically -relatively separated: [i] Renormalization 
group evolution (RGE), {ii) physical spectrum calculations (PS), {iii) effective potential 

^•^The program can be down-loaded from the node: http://lpm.univ-montp2.fr:7082/ djouadi/gdr.html 
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calculation with implementation of electroweak symmetry breaking (EWSB). The overall 
algorithm then reads as follows: choice of a model assumption/option choice of initial 
scale Qin (driven from input file suspectl2.in or from user's main code) —>■ RG evolution 
—>■ consistency of EWSB which involves the effective potential at one-loop (iterating until 
stability is reached) —>■ physical spectrum calculation: gauginos, sfermions, Higgses — >■ final 
masses and results (warning + comments as well) collected in file suspect. out. 

An important aspect of SUSPECT is a special attention given to the consistency of EWSB, 
which makes that not all of the scalar sector parameters are independent. [For the moment 
only the simplest constraints diVcs) / dvu,d = are included; the constraints from the absence 
of Charge and Color Breaking (CCB) minima will be implemented in a later version]. In 
particular, this is used to define different set of input/output scalar parameters. Although 
this resulting flexibility in the choice of input parameters is welcome, its actual implemen- 
tation is quite non trivial, which is payed by a slower CPU time. Moreover, one should keep 
in mind that it is often a main source of possible discrepancies with other similar task codes 
which implement EWSB in a different way. 

Another important ingredient of SUSPECT is the implementation of RG evolution, in 
different (loop) approximations. The RGE can be implemented (or not) by using different 
ichoice(l) input parameters. For instance, for ichoice(l)=0 one has the unconstrained 
MSSM with no RGE, i.e. the relevant input parameter are assumed to be at LOW scale. For 
ichoice(l) = 1, RGE in the unconstrained MSSM with non-universality and the inputs 
are assumed at high scale, except tan/3 to be given at low energies, ichoice(l) = 2: 
unconstrained MSSM with RGE bottom-up; the relevant input is set similarly as with 
ichoice(l)= 0, but the final output consists of all the soft parameters at the high scale. 
ichoice(l) = 10: minimal SUGRA model. 

For interfacing SUSPECTl . 2 with your main code, all the user has to control is the way to 
dialog between her/his "main" routine/program and the SUSPECTl. 2 subroutine, together 
with the precise meaning of the different "dialog" parameters, which are of two kinds: 

- The "physical" parameters, are those parameters that are either necessary input for a 
given model and/or running option, or the desired output. All such parameters are passed 
from the calling code to SUSPECT and back via specific COMMONS. By "physical" we mean 
either truly physical parameters such as masses etc [and that are generally the output of 
SUSPECT calculation], or MSSM basic parameters such as the SUSY and soft-SUSY breaking 
terms of the MSSM Lagrangian, that are generally input for the SUSPECT calculation. 

- The "control" parameters, whose different purpose is to choose various running options. 
There are three main "control" parameters appearing as arguments of the SUSPECT calling 
command: (i) iknowl sets some degree of control on various parts of the algorithm [=0 blind 
use, i.e. no control on any "algorithmic parameter, =1 more educated use, (ii) input setting 
control [=0 relevant parameters are read form suspect 12. in and =1 define the relevant 
inputs from your calling program] and [iii) ichoice for the choice of model parameters with 
ichoice(l) discussed above for the RGE and ichoice (6) for the scalar sector input [=0 for 
fi, Ma inputs and =1 for Mjj^, M^^ as inputs]. 

All details on the main core SUSPECT routines, input and output parameters as well as 
physical and control parameters can be found on the web site and in Ref. [Q. 
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5 SUSYGEN 



SUSYGEN2 [|3 is a Monte Carlo event generator for the production and decay of supersym- 
metric particles and has been initially designed for e~^e~ colliders. It has been extensively 
used by all four LEP experiments to simulate the expected signals. It includes pair produc- 
tion of charginos and neutralinos, scalar leptons and quarks. It offers also a possibility to 
study the production of a gravitino plus a neutralino within GMSB models and the produc- 
tion of single gauginos if one assumes R-Parity to be broken. 

All important decay modes of SUSY particles relevant to LEP energies have been im- 
plemented, including cascades, radiative decays and R-Parity violating decays to standard 
model particles. The decay is included through the exact matrix elements. The lightest su- 
persymmetric particle (LSP) can either be the neutralino x?? the sneutrino z> or the gravitino 
G in R-parity conserving models, or any SUSY particle if R-parity is violated. 

The initial state radiative corrections take account of pr/pL effects in the Structure Func- 
tion formalism. QED final state radiation is implemented using the PHOTOS fl^ library. An 
optimized hadronization interface to JETSET is provided, which also takes into account 
lifetimes of sparticles. Finally, a widely used feature of SUSYGEN2 is the possibility to perform 
automatic scans on the parameter space through user friendly ntuples. 

Recently SUSYGEN2 has been upgraded to SUSYGEN3 0] in order to adapt to the needs of 
the next generation of linear colliders, but also in order to extend its potential to supersym- 
metric particles searches at e~p colliders (e.g HERA) and hadronic colliders (e.g Tevatron or 
LHC) . The main new features relevant for linear colliders are the inclusion of beamstrahlung 
through an interface to CIRCEE |T^, the full spin correlation in initial and final states, the 
inclusion of CP violating phases and the possibility to have an elaborate calculation of the 
MSUGRA spectrum through an interface to SUSPECT |[. 

a) Mass spectrum calculation: 

SUSYGEN2 offers different frameworks for the mass spectrum calculation. One can first 
assume the different mass parameters entering in the MSSM: Mi, M2 and M3 the gaugino 
mass parameters, /i, the Higgsino mass mixing parameter, the scalar fermions masses Mj^ 
and Mj:^, the trilinear mixing parameters At and Aj. to be free. This gives the so called 
"unconstrained MSSM" . Another approach to the mass spectrum calculation is based on the 
supergravity inspired models. In this case the soft breaking mass parameters are assumed 
to be universal at the GUT scale reducing the number of parameters to mi/2, the common 
gaugino mass parameter, mo, the common sfermion mass parameter, the sign of fi, tan/5, 
the ratio of the two vacuum expectation values of the two Higgs doublets, Aq, the common 
trilinear couplings. All these parameters are defined at the GUT scale. 

In SUSYGEN3, one can keep the approach used in SUSYGEN2. In this case, only mo is 
defined at the GUT scale and the sfermion masses are evolved from the GUT scale to the 
electroweak (EW) scale according to the formulae given in appendix of Ref . . The other 



parameters Mi, fx, At, A^ and Aj- are defined at the EW scale and mixing of the third 
generation sfermion is taken into account through the parameters At, A^ and A^-. SUSYGEN3 
offers now the possibility to do a better treatment of the mass spectrum calculation within 
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mSUGRA through an interface to the SUSPECT program . In practice, if the flag SUSPECT 
is set to TRUE in the input data card which fixes the model, the entire mass parameters at the 
EW scale will be derived from these at the GUT scale (mi/2, mo, sign of /i, Aq and tan/?). 

b) Beam polarization and spin correlations 

Since one expects high luminosities for the next generation of linear colliders (e.g. ~ 500 
fb^^ for the TESLA project), one can use beam polarization to reduce the standard model 
backgrounds and use the polarization dependence of the cross sections to study specific SUSY 



parameters. Moreover, as it has been stressed by several authors [21|, spin correlations play a 
major role in the kinematic distributions of final particles. To fulfill these two requirements, 
the "helicity amplitude method" was used for the calculation of the different Feynman 
amplitudes for production and decay, in order to obtain full spin correlation. Since such 
amplitudes involve products and contractions of fermionic currents, two basic functions, 
namely the B and Z functions were defined through: 

^Ai,A2(Pl>P2) = Ux,{pi,mi)PxUx2{p2,m2) (1) 
^a^'a2,A3,A4(P1'P2,P3,P4) = [ux,{pi,mi)-f^'PxUx2{p2,m2)] [uxAP3,mshf,Px'Ux,{p4,m4)] 

where Px stands for one of the two chiral projectors Pi or Pr and ux{p, m) denotes the posi- 
tive energy spinor solution of the Dirac equation for a particle of helicity A, four momentum 
p and mass m. The decomposition of the bispinors ux{p, ttl) in terms of the massless helicity 
eigenstates uJx{k) yields simple analytical expressions for the B and Z functions. The ampli- 
tude is then factorized in terms of these basic building blocks; this fact permits compact and 
transparent coding and speed of calculation. The masses are not neglected in any stage of 
the calculation. For gaugino productions and decay, we use the "widthless approximation" . 
For instance, the calculation of the cross section associated to e"'"e~ — > X2X1 ~^ XiXi^~^^~ is 
done as follows: the total amplitude associated to a given helicity configuration of the differ- 
ent particles is approximated by the product of the amplitude associated to the production 
of the two neutralinos A^(e+e~ — > X2X1) with the amplitude corresponding to the decay of 
the next to lightest neutralino A4{x2 ~^ x\^^^~)- The remnant of the propagator squared 
of X2 is approximated by a factor given by 87r''/(m^or^o). The phase space integration is 
done through the multichannel method 



c) Including phases in SUSY searches: 



In the MSSM, there are new potential sources of CP non-conservation |2^. Complex 
CP violating phases can arise from several parameters present in the MSSM Lagrangian: 
the higgs mixing mass parameter /i, the gaugino masses Mj, the trilinear couplings Ai. 
Experimental constraints on these CP violating phases come from the electric dipole moment 
of the electron and the neutron. Since in SUSYGEN3 all the couplings, the different mass 
parameters /x. Mi, and the trilinear couplings A^, At and A^ have been assumed to be 
complex by default [E^, the introduction of phases in the gaugino and sfermion sector for 



masses as well for cross sections has been straightforward. 
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6 CompHEP 



CompHEP[^ p9| is a package for automatic calculations of decay and production processes in 
the tree-level approximation in the framework of arbitrary gauge models of particle interac- 
tions. The main idea prescribed into CompHEP, is to make available passing on from the basic 
Lagrangian to the final distributions efficiently with a high level of automation. CompHEP is 
a menu-driven system. The codes and the manual are available on the web site: 
http://theory.npi.msu.su/~comphep (mirror on http://www.ifh.de/~pukhov). 

The present version has four built-in physical models. Two of them are the Standard 
Model in the unitary and 't Hooft-Feynman gauges. The user can change particle inter- 
actions and model parameters and introduce new vertices, thus creating new models. Fur- 



thermore, in the framework of the CompHEP project, a program LanHEP was created to 
generate CompHEP model files as will be discussed below. 

The CompHEP package consists of two parts, a symbolic and a numerical one. The symbolic 
part is written in the C programming language and produces FORTRAN and C codes for squared 
matrix elements which are used in the numerical calculation later on. There are two versions 
of the numerical part, a FORTRAN and a C one, with almost equal facilities. The C version 
has a more comfortable interface but it does not possess an option to generate events and 
does not perform calculations with quadruple precision. 

The symbolic part of CompHEP allows the user to: 

- Select a process by specifying incoming and outgoing particles for the decays 1 — > n 
(< 6) and the production mechanisms 2 — n (< 5). 

- Generate Feynman diagrams, display them, and generate squared diagrams. 

- Calculate analytical expressions corresponding to squared diagrams, save them in 
REDUCE and MATHEMATICA forms for further symbolic manipulations. 

- Generate optimized FORTRAN and C codes for the squared matrix elements for further 
numerical calculations. 

The numerical part of CompHEP allows to: 

- Convolute the squared matrix element with structure functions (for proton and an- 
tiproton, electrons and photons). 

- Modify physical parameters (energy, charges, masses etc.) involved in the processes. 

- Select the scale for evaluation of as and parton structure functions. 

- Introduce various kinematical cuts. 

- Define the phase space parameterization and introduce a phase space mapping in order 
to smooth sharp peaks for effective Monte Carlo integration. 

- Perform Monte-Carlo integrations by VEGAS via the multichannel approach [^]. 

- Generate events and make distributions with graphical and LaTeX outputs. 

In the QCD part of these Proceedings, one can find more details on CompHEP options, in 
particular the handling of the QCD aspects and the discussion of the automatic computation 
of processes with multiparticle final states. The CompHEP package has been used in several 
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studies performed at this Workshop, in particular in the Higgs working group. Examples 
are: Higgs boson searches in the 77+jet channel at the LHC [Sec. 2] and generation of events 
for associated production of light stops with Higgs bosons [Sec. 4] . 

During this Workshop, a new algorithm was proposed for the treatment of the first and 
second generation quarks through the single generation of generalized "up" and "down" 
quarks [^. This algorithm neglects the masses of these quarks and their mixing with third 
generation quarks. It is based on a rotation of the S-matrix in flavor space and move the 
CKM matrix elements from diagrams to distribution functions. The complete set of new 
rules was derived for a correct counting of the convolution with different parton distributions 
for quarks of the first/second generations. Each rule corresponds to a gauge invariant subset 
of diagrams; see also ||3^. This technique allows to reduce significantly the number of sub- 
processes contributing to the same physical final state, especially for hadron colliders. It was 
realized in the CompHEP version installed at CERN (/af s/cern.ch/cms/physics/COMPHEP). 

Developments were also made during this Workshop for the implementation of SUSY 
models in CompHEP; some of them concern the Higgs sector. To derive the MSSM description 



for CompHEP one can use the LanHEP [^ program which allows to generate the Feynman 
rules from the Lagrangian input in compact forms close to the ones given in textbooks [e.g. 
Lagrangian terms can be written with summation over indices and using compact expressions 
such as covariant derivatives and strength tensors for gauge fields]. There are given in 
terms of two-component spinors and with the superpotential formalism. The output for the 
Feynman rules is in LaTeX format and in the form of CompHEP model files. For the MSSM 
Lagrangian, the complete description given in Ref. [^ is used, together with two extensions: 
vertices with R-parity violation and the light gravitino scenario in GMSB models. 

It is known that Higgs boson masses in the MSSM are significantly affected by radia- 
tive corrections. To compute these corrections, the two-Higgs doublet model potential p5| 
technique is exploited. This potential is parametrized by 7 variables, A1...A7, for which an- 
alytical formulae given by in M. Carena et al. in Ref. [Q] are implemented. CompHEP allows 
to calculate arbitrary processes within the given physical model. Thus, one has to deal 
with the Xi variables rather than with the set of Higgs boson masses only. However, one 
can set the Higgs boson masses as input parameters, but the Aj are derived after and the 
model is changed correspondingly preserving gauge invariance. An interface is made with the 
FeynHiggs program [|1^] [used as an external library], thus providing an option to evaluate 
the CP-even Higgs boson masses in the most up-to-date way. 

The number of independent parameters in the MSSM can be reduced by the imple- 
mentation of the mSUGRA or GMSB models. More specifically, the soft SUSY-breaking 
parameters, gaugino and sfermion masses as well as trilinear couplings, are computed from 
the input parameters. It is possible to use the ISASUSY package [^ for the calculation of 
these soft SUSY-breaking parameters [as well as the CP-odd Higgs boson mass; the CP- 
even Higgs masses can be calculated by FeynHiggs]. The masses of the sparticles are then 
calculated by CompHEP from the formulae used in the unconstrained MSSM. SUSY models 
for CompHEP with the FeynHiggs and ISASUSY options included, can be obtained from the 
web site: http://theory.npi.msu.su/~semenov/mssm.html 
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